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CN ; Abstract 

I ' We show that the QCD vacuum (without dynamical quarks) is a dual supercon- 

Q\ • ductor at least in the low-energy region in the sense that monopole condensation does 

^ . really occur. In fact, we derive the dual Ginzburg-Landau theory (i.e., dual Abelian 

O ■ Higgs model) directly from the SU(2) Yang-Mills theory by adopting the maximal 

Abelian gauge. The dual superconductor can be on the border between type II and 
type I, excluding the London limit. The masses of the dual Abelian gauge field is 
expressed by the Yang-Mills gauge coupling constant and the mass of the off-diagonal 
^ . gluon of the original Yang-Mills theory. Moreover, we can rewrite the Yang- Mills the- 

ory into an theory written in terms of the Abelian magnetic monopole alone at least 
'k/. \ in the low-energy region. Magnetic monopole condensation originates in the non-zero 

^ I mass of off-diagonal gluons. Finally, we derive the confining string theory describing 

the low-energy Gluodynamics. Then the area law of the large Wilson loop is an im- 
mediate consequence of these constructions. Three low-energy effective theories give 
the same string tension. 
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1 Introduction 

The main aim of this paper is to discuss how the dual superconductor picture for 
explaining quark confinement is derived directly from Quantum Chromodynamics 
(QCD). It is believed that quark and gluon are confined into the inside of hadrons by 
the strong interaction described by QCD and that they can not be observed in isola- 
tion. Quark confinement can be understood based on an idea of the electro-magnetic 
duality of ordinary superconductivity, i.e., the dual superconductor picture|[l|]. In this 
picture, the color electric fiux can be excluded from QCD vacuum as a dual super- 
conductor (the dual Meissner effect), just as the magnetic field can not penetrate 
into the superconductor (the so-called Meissner effect). In this context, the dual is 
used as implying electro-magnetic duality. However, QCD is a non-Abelian gauge 
theory and its gluonic part is described by the Yang-Mills theory, whereas the usual 
superconductivity is described by the Abelian gauge theory. Therefore we must make 
clear the precise meaning of dual superconductivity in QCD. It is possible to assume 
that the diagonal component of color electric field can be identified with the dual 
of magnetic field in the ordinary superconductivity. If so, we must answer what is 
the role played by the off-diagonal component of gluons? Finally, we must show that 
QCD vacuum has a tendency to exclude the (diagonal component of) color electric 
field. If a pair of heavy quark and anti-quark is immersed in the QCD vacuum, the 
color electric fiux connecting a pair is squeezed into the tube-like region, leading to 



the formation of QCD (gluon) string between a quark and an anti-quark. As a re- 
sult, the interquark potential V{R) is proportional to the interquark distance R. To 
separate a quark from an anti-quark, we need infinite energy. In this sense, quark 
confinement is achieved. 

The dual superconductor picture of QCD vacuum is based on the following as- 
sumptions: 

1. Existence of magnetic monopole: QCD has magnetic monopole. 

2. Monopole condensation: Magnetic monopole is condensed in QCD. 

3. Infrared Abelian dominance: Charged gluon (i.e., off-diagonal gluon) can be 
neglected at least in the low-energy region of QCD. 

4. Absence of quantum effect: Classical configuration, e.g., the magnetic monopole, 
is dominant. Hence quantum effect can be neglected. 

Each of these statements should be explained based on QCD to really confirm the dual 
superconductor picture of QCD vacuum. As a first step, a prescription of extracting 
the field configuration corresponding to magnetic monopole was proposed by 't Hooft 
0. This idea is called the Abelian projection. Immediately after the proposal of this 
idea, infrared Abelian dominance was suggested by Ezawa and Iwazaki [Q]. However, 
recent simulations revealed that the Abelian dominance is not necessarily realized in 
all the Abelian gauges, see the review [^]. To author's knowledge, the best covariant 
gauge fixing condition realizing Abelian dominance is given by the maximal Abelian 
(MA) gauge proposed by Kronfeld et al. f^. In fact, the infrared Abelian dominance 
is first confirmed in MA gauge based on Monte Carlo simulations by Suzuki and Yot- 
suyanagi 0. Subsequent simulations have also confirmed the monopole dominance 
in low-energy QCD[|^. An analytical derivation of the dual Ginzburg-Landau (DGL) 
theory was tried by Suzuki [|| by way of Zwanziger formulation by neglecting the 
off-diagonal components of gluon fields by virtue of the Abelian dominance and by 
assuming condensation of magnetic monopoles. However, the DGL theory is not yet 
derived directly from the underlying theory, i.e. QCD, since the Abelian dominance 
and monopole condensation themselves must be derived in the same framework of 



the theory. According to the recent Monte Carlo simulations |T3] , the presumed dual 
Ginzburg-Landau theory is on the border of type II and type I. 

In this paper we discuss how to derive the DGL theory as a low-energy effective 
theory (LEET) of QCD. The LEET is not unique and we can derive various LEET's. 
Of course, they should be equivalent to each other, if they are to be derived directly 
from QCD. Even if a specific LEET of QCD is assumed, the parameters included in 
the LEET can be adjusted so as to meet the data of experiments or Monte Carlo 



simulation on a lattice. In the first paper |TT| of a series of papers []T2|, [13], |T^, |T3|, |K 
on the quark confinement in Yang-Mills theory in the Abelian gauge, the author has 
given a scenario of deriving the dual superconductivity in low-energy region of QCD 
and demonstrated that the DGL theory, i.e., the dual Abelian Higgs (DAH) theory 
in the London limit can be derived from QCD, provided that the condensation of 



magnetic monopole takes place. In this scenario, the mass rrib of the dual gluon (dual 
Abelian gauge field) is given directly by the non-zero condensate of the magnetic 
monopole current /c^, i.e., {k^k^) ^ 0. Since the DGL theory is a LEET written 
in term of only the Abelian (diagonal) component extracted from the original non- 
Abelian field, the treatment of the off-diagonal component is crucial for deriving the 
LEET and also for giving reasonable interpretation of the result. In the previous paper 



111 , all the off-diagonal components are integrated out to write down the LEET in 
terms of the diagonal components alone. The resulting LEET is an Abelian gauge 
theory preserving a characteristic feature of non-Abelian gauge theory in the sense 
that the the j3 function for the coupling constant has the same form as the original 



Yang-Mills theory, exhibiting the asymptotic freedom [|1^, ^. For this procedure to 
be meaningful, the mass Ma of the off-diagonal components must be heavier than 
the mass of the diagonal components. This procedure can be justified based on the 
Wilsonian renormalization group or the decoupling theorem [^]. As a result, the 
LEET is valid in the energy region below M^. 

The aim of this paper is to show that the dual superconductivity can be derived 
at least in the low-energy region of Gluodynamics (i.e., the gluonic sector of QCD). 
The most important ingredient in deriving the dual superconductivity is the existence 



of non-zero mass of the off-diagonal gluons. Recently, it has been shown [19, 2C] that 
the off-diagonal gluons and off-diagonal ghosts (anti-ghosts) acquire non-zero masses 
in Yang-Mills theory in the modified MA gauge [^]. This result strongly supports 
the Abelian dominance in low-energy Gluodynamics. The modified MA gauge was 



already proposed by the author from a different viewpoint in the paper []12[. A 
remarkable difference of MA gauge from the usual Lorentz type gauge lies in the 
fact that MA gauge is a nonlinear gauge. In order to preserve the renormalizabil- 
ity of Yang-Mills theory in the MA gauge, it is indispensable to introduce quartic 
ghost interaction term as a piece of gauge fixing term [^. The modified MA gauge 
fixes the strength of the quartic ghost interaction by imposing the symmetry, i.e., 
orthosymplectic group 05'p(4|2). The implications of the 0Sp{A\2) symmetry have 
been discussed in the previous papers [0, |TB|, |T^. The attractive quartic ghost in- 
teraction causes the ghost-anti-ghost condensation. Consequently, the off-diagonal 
gluons become massive,^ whereas the diagonal gluons remain massless in this gauge. 
Since the classical Yang-Mills theory is a scale invariant theory, the mass scale must 
be generated due to quantum effect. In this paper, we show that the monopole con- 
densation does really occur due to existence of non-zero mass of off- diagonal gluons. 
The off-diagonal gluon mass also provides the mass of dual gauge field in the DGL 
theory. 

It should be remarked that in the conventional approaches the off-diagonal com- 
ponents are completely neglected from the beginning in deriving the effective Abelian 
gauge theory by virtue of the Abelian dominance. As a result, the conventional ap- 
proach can not predict the physical quantities without some fitting of the parameters 
introduced by hand. The purpose of this paper is to bridge between the perturbative 



^The off-diagonal gluon mass Ma in the MA gauge has been calculated on a lattice by Amemiya 
and Suganuma ||, Ma = 1.2GeV for SU{i). 



QCD in the high-energy region and the low-energy effective Abehan gauge theory. 
Consequently, the undetermined parameters in the low-energy effective theory can in 
principle be expressed by the parameters of the original Yang-Mills theory, i.e., the 
gauge coupling constant g and the renormalization group (RG) invariant scale Aqcd- 
In this paper we pay attention to the vacuum expectation value (VEV) of the 
Wilson loop operator iy(7[^], which is written in the framework of the functional 
integration as 

{Wc[A])ym = Zyl, J dfiexpi^t J d^xC'fM]Wc[A], (1.1) 

where CyM is the total Yang-Mills Lagrangian, i.e., Yang-Mills Lagrangian Cym 
plus the gauge fixing (GF) term Cgf+fp including the Faddeev-Popov (FP) ghost 
term in the modified MA gauge. Our derivation is based on the Becchi-Rouet-Stora- 
Tyutin [^(BRST) formulation of the gauge theory. Hence, d^ is the integration mea- 
sure which is ERST invariant, d/iyM '■= 'DA^'DBVC'DC , where B is the Nakanishi- 
Lautrap (NL) Lagrangian multiplier field. We regard the Wilson loop as a special 
choice of the source term in general Yang-Mills theory. The Wilson loop is used as 
a probe to see the QCD vacuum. The shape of the Wilson loop is arbitrary at this 
stage. 

We seek other theories (with the corresponding source term) which is equivalent 
to the original Yang-Mills theory at least in the low-energy region in the sense that it 
gives the same VEV of the Wilson loop operator as the original Yang-Mills theory for 
large loop C. From this viewpoint, we obtain three LEET's, i.e., DGL theory (DAH 
theory), magnetic monopole theory and confining string theory. For this derivation 
to be successful, the existence of quantum corrections coming from off-diagonal com- 
ponents of gluons and ghost (and anti-ghost) is indispensable. Without quantum 
corrections, we can not derive magnetic monopole condensation. This is plausible, 
since it is the quantum correction that can introduce the scale into the gauge theory 
which is scale invariant at the classical level. 

This paper is organized as follows. In section ^, we give a strategy (steps) of 
deriving the LEFT of Gluodynamics. In fact, we demonstrate that at least the 
London limit of the DGL theory (i.e., DAH model) of type II can be obtained as a 
very special limit of the resulting LEFT of Gluodynamics. However, our method is 
able to treat more general situation, not restricted to the London limit. Rather, our 
result suggests that the DGL theory can not be in the London limit. 

In section ^, a more general case of the DGL theory is discussed. It is shown that 
the LEFT of the supposed DGL theory agrees with the LEFT derived from the Yang- 
Mills theory according to the above steps, in the energy region less than the mass 
iTiH of the dual Higgs mass. This way of showing equivalence is a little bit indirect. 
The reason is as follows. We know that the DAH model or the DGL theory is a 
renormalizable theory (within perturbation theory in the magnetic coupling constant 
gm = 4:T^/g) and hence is a meaningful theory at arbitrary energy scale. On the other 
hand, we know that the high energy region of Gluodynamics is correctly described 
by the non-Abelian Yang-Mills theory with asymptotic freedom. Therefore, the DGL 
theory is at best meaningful only in the low-energy region in this context, although it 



is a renormalizable theory. In this sense, we must be careful in saying that the DGL 
theory is regarded as a LEET of Gluodynamics. 

In section ^ we try to estimate the neglected terms in the derivation of the LEET 
based on the large N argument and the decoupling theorem. It is shown that in the 
large N expansion the higher-order terms are suppressed by iV~^ compared to the 
leading order term. Thus the LEET obtained above is considered as the leading-order 
result of the large A^ expansion. 

The LEET of Gluodynamics is not unique. In fact, we can obtain various LEET's 
which are equivalent to each other. Other LEET's can be more convenient for the 
purpose of calculating some kinds of physical quantities. 

In section ^ we derive the magnetic monopole action as another LEET. The mag- 
netic monopole action SMp[k] is written entirely in terms of the magnetic monopole 
current k. 



"p.- 



{Wc[A])ym = Z]^lpjvk^exp i^t J d^x{CMp[k] + fc^S'^)} , (1.2) 

by introducing the four- dimensional solid angle S^. We show that the VEV of the 
Wilson loop exhibits area law decay for large loop C, 

{Wc[A])ym = exp{-astA{C)}, (1.3) 

where A{C) is the area of the minimal surface S bounded by C. The string tension 
ast is explicitly obtained and it agrees with that predicted by the DGL theory.. Then 
we show that the magnetic monopole condensation really occurs in the low-energy 
region in the sense that the current- current correlation at the same space-time point, 
(k^kfj), has non- vanishing expectation value, 

{k^k>')YM ■■= Zllp J Vk^ exp |z j d^xCMp[k]^ k^k^ ^ 0. (1.4) 



The result is consistent with that of the monopole action on a lattice p^. 

In section ^ we derive a string theory which is equivalent to the LEET's derived 
above, 

{Wc[A])ym = Z;,'Jvx^{a)exp{iS,,[x]}. (1.5) 

The action Scs [x] of the confining string is equal to the Nambu-Goto action for the 
world sheet 5* with the Wilson loop C as the boundary. The string tension o"cs in the 
Nambu-Goto action is the same as ast evaluated by the monopole action derived in 
section ^. The Wilson loop exhibits area law with the string tension agt- Therefore, 
the obtained string theory is regarded as a low-energy limit of the so-called confining 



string theory proposed by Polyakov p5| 



In section ^ we discuss what values of the parameters in the LEET's should be 
chosen to reproduce the numerical results. 

The final section ^ is devoted to summarizing the result obtained in this paper 
and discussing the future directions of our investigations. The details of calculations 
are collected in Appendices, together with useful formulae. 



2 Dual superconductivity in low-energy Gluody- 
namics 

2.1 Conventions 

The gauge potential Afj, is written as 

A,{x) := A^{x)T^ iA = l,---,N'-l), (2.1) 

where the generators T^{A = 1, ■ ■ ■ , A^^ — 1) of the Lie algebra Q of the gauge group 
G = SU{N) are taken to be Hermitian satisfying [T^, T^] = if^^'-^'T'^ and normalized 
as tr(T'^T'^) = ^S^^. For a closed loop C, we define the Wilson loop operator W{C) 

by 



W{C) = tr 



Pexp < ig j> dx^A^{x] 



Ml), (2.2) 



where P denotes the path-ordered product and g is the Yang- Mills coupling constant. 
We begin with the vacuum expectation value (VEV) of the Wilson loop operator 
W{C) in the Yang-Mills theory with a gauge group SU{N) defined by the functional 
integral, 

{W{C))ym := Zyl, J dfiYMexp{zSpl)W{C), (2.3) 

where Zym is a normalization factor (or a partition function) to guarantee {1)ym = 1 
and Sym is the total action obtained by adding the gauge fixing (GF) and Faddeev- 
Popov (FP) ghost term Sgf+fp to the Yang-Mills action Sym, 

^YM = ^YM + Sgf+fp- (2.4) 

The Yang- Mills action is of the usual form, 

Sym = j d\^:F^,{x):F^'^\x), (2.5) 

where J^^y{x) is the field strength defined by 

^^.u{x) := T^,{x)T^:=d^A,{x)-d,A^{x)-ig[A^{x),A,{x)]. (2.6) 

The GF+FP action Sgf+fp is specified later (see section pT^ ). Finally, dfiYM is the 
integration measure, 

d/jYM := VA'^VB^VC^VC^, (2.7) 

which is invariant under the Becchi-Rouet-Stora-Tyutin (BRST) transformation, 

5bA^.{x) = V^[A\C{x) := df,C{x) - ^g\A^,[x),C{x)\, 
bBC{x) = i-[C{x),C{x)], 

6bC{x) = iB{x), 

5bB{x) = 0, (2.8) 

where B is the Nakanishi-Lautrap (NL) field and C (C) is the ghost (anti-ghost) field. 

6 



2.2 Step 1: Non-Abelian Stokes theorem for the Wilson loop 



We make use of a version of the non-Abelian Stokes theorem (NAST) |2^, |T^, |2^ to 
rewrite the Wilson loop operator in terms of the diagonal components. This version 
of NAST was first obtained by Diakonov and Petrov for SU{2) |]2B|. It is possible to 
generalize their result to SU{N) {N > 3). 

Theorem: ||27]| For a closed loop C , we define the non-Abelian Wilson loop operator 
by 



Wc[A] = trjPexp ig d) i 



LLtLi '-'^'LiK'^ } 



/tr(l) 



where P is the path- ordered product. Then it is rewritten as 



(2.9) 



Wc[A] = J dficiO exp ig j 



ig f dx^a^^{x) 



c^Atc(Oexp 



'i9J/a^'fi. 



x] 



, (2.10) 



where 




(2.11) 



anc 



I flU 



(x) ■.= d^al{x) -d^al{x). 



(2.12) 



Here |A) is the highest-weight state of the representation defining the Wilson loop and 
the measure dficiO ^■^ the product measure along the loop C , d^ciC) = Hx^c ^f^i^i^)) > 
where dfi{C,) is the invariant Haar measure on G/H with the maximal stability group 
H . The maximal stability group H is the subgroup leaving the highest-weight state 
invariant (up to a phase factor), i.e., 



g\K) = ih\K) = i\A)e 



i(l>(h) 



(2.13) 



for ^ G G/H and h & H. It depends on the group G and the representation in 
question. 

For G = SU{2), the H is given by the maximal torus subgroup H = f/(l) ir- 
respective of the representation. Hence G/H = SU{2)/U{1) = GP^ = Fi. For 
G = SU{N){N > 3), however, H does not necessarily agree with the maximal torus 
group H = U{1)^~^ depending on the representation. For G = SU{3), all the 
representations can be classified by the Dynkin index [m, n]. If m = or ra = 0, 
H = U{2) and G/H = GP"^. On the other hand, when m 7^ and n 7^ 0, 
H = f/(l)2 = [/(I) X U{1) and G/H = F2. Here GP'' is the complex projective 
space and Fn the flag space. This NAST is obtained by making use of the generalized 



coherent state. For details, see the reference [27 



Note that ffj,i, is not equal to the diagonal component of JF^jy. 



For the fundamental representation, the expression ( p.ll| ) is greatly simplified as 



(A|(.--)|A) 
Therefore, 



2tr[7^(--.)], 



Ti = -diag 



A^-1 -1 



-1 

aT 



2^ 



for G = 

4 + -^K 



SU{2) 

for G = SU{3). 



(2.14) 

(2.15) 
(2.16) 



This implies that the non-Abelian Wilson loop can be expressed by the diagonal 
( Abelian) components. This is suggestive of the Abelian dominance in the expectation 
value of the Wilson loop. 

The monopole dominance in the Wilson loop is also expected to hold as shown 
follows. We can rewrite /^j, in the NAST as 



f^x) = d,[n''ix)ACix)] - d4n^ix)A^ix)] - -r^'^n^{x)d,n^{x)d,n^ {x), (2.17) 
where 

n^{x)T'^ = i{x)ni\x) = u{x)nu\x). (2.18) 

The /^^ is invariant under the full G gauge transformation as well as the residual H 
gauge transformation. (Indeed, we can write a manifestly gauge invariant form, see 
ref. ||2^ .) The /^j, is a generalization of the 't Hooft-Polyakov tensor of the non-Abelian 
magnetic monopole, if we identify n"^[x) with the unit vector of the elementary Higgs 
scalar field in the gauge-Higgs theory: 



n^ix) 



^\x):=<p\x)/\<P{x)\. 



(2.19) 



This implies that n'^[x) is identified with the composite scalar field and plays the 
same role as the scalar field in the gauge-Higgs model, even though QCD has no 
elementary scalar field. This fact could explain why the QCD vacuum can be dual 
superconductor due to magnetic monopole condensation. 

By introducing the magnetic monopole current k hj k := 6 * f (see Appendix ^, 
we have another expression. 



(2.20) 



where A is the Laplacian A := d6 + 6d and T is a two-form determined by the surface 
element dS of the surface spanned by the Wilson loop C. The derivation is given in 
Appendix |B[ Hence, the Wilson loop can also be expressed by the magnetic monopole 
current fc^. The above results hold irrespective of which gauge theory we consider. 
In the case of SU{2), the Wilson loop in an arbitrary representation is written in 




the form [26, 



Wc[A] 



dficiO^W 



igj j> dx^a^^ 



X 



(2.21) 



where a^ is the Abelian gauge field (or the diagonal component of A^^) defined by 

tr{a3[e(x)^^(x)e^(x) + ^g''ax)^,e{x)]}, (2.22) 



4(x) 



for an element ^ G G/H = SU{2)/U{1) ^ S^ = CP^. Here J is a character which 
distinguishes the different representation defining the Wilson loop, J = i,l,|,.... 
Moreover, the use of the usual Stokes theorem leads to 



(2.23) 




where /L(a^) is the Abelian field strength defined by /^ 



fiu 



X 



d^al{x) 



d^a^^ix) 



and Sc is an arbitrary two-dimensional surface with the loop C as the boundary. 
Here it should be remarked that ff;_^{x) is invariant under the full G = SU{2) gauge 
transformation as well as the residual H = U{1) gauge transformation, since it has the 
same form as the usual 'tHooft-Polyakov tensor describing the magnetic monopole, 

see 0. 



2.3 Step 2: Cumulant expansion 

Now we specify the gauge theory in terms of which the VEV of the Wilson loop 
operator is evaluated. We consider the Yang-Mills theory with gauge fixing term. 
The gauge fixing is discussed in the next step. 

By applying the cumulant expansion to the VEV of the Wilson loop. 



{W{C))ym= /rf/ic(0(exp 



tgj dS^^r^ 



Sc 



fiu 



X 



YM 



(2.24) 



the VEV of the exponential is replaced by the exponential of the connected expecta- 
tion as 



{WiG))YM= /rf/ic(Oexp 



2 

\.P f dS^^{x)f dSP'^{y){fl^{x)flM)yM + -- 

^ J Sc ^ Sc 

(2.25) 
where we have used {f^ui^))YM = and ■ ■ ■ denotes the higher-order cumulants. 

The cumulant expansion is a well-known technique in statistical mechanics and 
quantum field theory. In what follows, we will neglect the higher-order cumulants 
in (|2.25|) as in the analysis of the stochastic vacuum model.0 The validity of this 
approximation, i.e, the truncation of the cumulant expansion, can be examined by 
Monte Carlo simulations on a lattice, as performed for the stochastic vacuum model, 
see [0. In the framework of our approach, this approximation can be justified in 
the sense that the approximation is self-consistent within the APEGT derived below. 



See section KA\ for more details. 



•^In the non-perturbative study of QCD, the cumulant expansion is extensively utilized by the 
stochastic vacuum model (SVM) ^^ where the different version of the non- Abelian Stokes theorem 
is adopted. In the SVM, the approximation of neglecting higher order cumulants is called the 
bilocal approximation. The validity of bilocal approximation in SVM was confirmed by Monte Carlo 
simulation on a lattice |B(j, aw. The author would like to thank Dmitri Antonov for this information. 



2.4 Step 3: Maximal Abelian gauge fixing 



First of all, we define the decomposition of the gauge potential into the diagonal and 
off-diagonal components, 



A,{x) = A^ix)T^ = aUx)T + Al{x)T\ 



(2.26) 



where T^ & H and T" G ^ — 7Y with H being the Cartan subalgebra. As a gauge 
fixing condition for the off-diagonal component, we adopt the modified version of the 
maximal Abelian (MA) gauge proposed by the author [jl2|. 



S, 



GF+FP 



I 



ct^xiSs^B 



\aI{x)A^\x) - '^tC\x)C\x) 



(2.27) 



where a corresponds to the gauge fixing parameter for the off-diagonal components, 
since the explicit calculation of the anti-BRST transformation 5b yields 



Sgf+fp = - d xi6B 



fy \ ry — — ry — — 

" ' ■ "abi/^a/^h/^i ■ ^ rabc/^a/^b/^c 



C \ D^[a]A'' + -3^ - i-gT'C'C'C' - i-gr^CC'C 



(2.28) 



In order to see the effect of ghost self-interaction, we take 
SnF-i-FP = — I d!^xi5 



>B 



C \ D^lajAf" + -bV - I'^gT^'C'C^C' - I'^gf^^C'C^C^ 



(2.29) 

where we must put C, = a to recover Eq.( p.27| ). The most general form of the MA 
gauge was obtained by Hata and Niigata [^ . 

By performing the BRST transformation explicitly, we obtain 

Sgf+fp = J d'x{B''D,[afA>'' + ^B^'B'^ 

+iC''D^[a]''%gr'^^A^'^C^) + iC''gr''\D^'[a]^'A''^)C' 

_i__5_ 2 rabe rcderia/^b/^cz-yd , 2. „2 rabc raid/^b/^c/^i/^d \ 2_„ fibc- r)b /~ia r^c 



8' 



4^ 



-CgP^HB^C^C' + 'kg'^pbijcdi^a^b^c^dy 



In particular, the SU(2) case is greatly simplified as 

Sgf+fp = J d''x{B''D^[afAf^' + ^B'^B'^ 

+iC''D^[a]'''D^'[aY^C^ - tgh'"^e^''C''C^A'''A 
+iC^ge^\D^[a]''''Al)C^ 

-Cge^HB^C^C^ + ^g^e''''e^'^C''C^C^C'^}. 



(2.30) 



(2.31) 
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Integrating out the NL field B"^ leads to 

Sgf+fp = ld''x{-^iD^[ar'A'^'f + {l-C/a)tC^ge''\D^[a]'''A';^)C' 

j^^g2^ab^cd^a^bfjcfjd^_ (2.32) 

The advantages of the modified MA gauge ( p.27| ) are as follows. 

1. Sgf+fp is BRST invariant, i.e., 6bSgf+fp = 0, due to nilpotency of the BRST 
transformation 6"^ = 0. 

2. Sgf+fp is anti-BRST invariant, i.e., SbSgf+fp = 0, due to nilpotency of the 
anti-BRST transformation 6^ = 0. 

3. Sgf+fp is supersymmetric, i.e., invariant under the 0Sp{A\2) rotation among 
the component fields in the supermultiplet {A^, C, C) defined on the superspace 
{Xfj,,6,6). The hidden supersymmetry causes the dimensional reduction in the 
sense of Parisi-Sourlas. Then the 4-dimensional GF+FP sector reduces to the 
2-dimensional coset (G/H) nonlinear sigma model. See ref.[|l^ for more details. 

4. Sgf+fp is invariant under the FP ghost conjugation, 

C^^±C^, C^^^C^, B^^~B^, B^^~B^, (A^-^A^). 

(2.33) 
Therefore, we can treat C and C on equal footing. In other words, the theory 
is totally symmetric under the exchange of C and C. 

5. The Yang-Mills theory in the modified MA gauge (with the total action Sym + 
Sgf+fp) is renormalizable. The naive MA gauge 



Sgf+fp = - J d'xz6B[c^ [D^[a]A'' + ^bJ 
= f d^x{B''D^[a]^^A^'^ + -B^B^ 



(2.34) 



+iC^DM'"'D>'[aY^C^ - tgh'"^e'''C^C^A'''Af^ 



ya „^ab( n \ „Abc 



+iC''ge''\D^[aX^A^^)C-' (2.35) 

spoils the renormalizability, since radiative corrections induce (even for a = 0) 
the four-ghost interaction, 

Zi^g^e^h^'^CC^C^C'^, Z4c = 4iV-^ In -^, (2.36) 

{A-ny /io 

owing to the existence of the vertex —ig'^e"''^e'^^C°'C^A^'^A'l, see eq.(2.52) and 



Appendix B of the paper |TT|. This is because the MA gauge is a nonlinear 
gauge. For the renormalizability of the Yang-Mills theory in the MA gauge, 
therefore, we need the four-ghost interaction from the beginning. In fact, the 
renormalizability of the Yang-Mills theory supplemented with the four-ghost 
interaction was proved to all orders in perturbation theory pT| . 
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In order to completely fix the gauge degrees of freedom, we add the GF+FP term for 
the diagonal component a* to ( p.28|) or ( p.29|) : 



S, 



GF+FP 



— d xi6 



B 



P 



CM9X + §5^ 



(2.37) 



where we have adopted the gauge fixing condition of the Lorentz type, d^a''^ = 0. The 
choice of the modified MA gauge is essential in deriving the off-diagonal gluon mass. 

2.5 Step 4: Dynamical mass generation for ofF-diagonal com- 
ponents 

It is shown that the four-ghost self-interaction is indispensable for the renormalizabil- 



ity of Yang-Mills theory in the MA gauge. Moreover, it has been shown [IS, 20] that 
the attractive four-ghost interaction in the modified MA gauge causes the ghost-anti- 
ghost condensation and that this condensation provides masses for the off-diagonal 
gluons and off-diagonal ghosts and anti-ghosts. The massive off-diagonal fields do not 
propagate in the long distance. Therefore, we can neglect off-diagonal components in 
the low-energy or long-distance region, except for the renormalization of the remain- 
ing diagonal fields. This result strongly supports the infrared Abelian dominance 
conjectured by Ezawa and Iwazaki Pj. 

The dynamical mass generation of the off-diagonal components is understood 
based on the argument of Coleman- Weinberg type. See the paper [^. The off- 
diagonal gluon propagator is given by 



mx)Al(y)) 



D»t(A:):=rD^.(A;), D,,{k) :- 



" '^ ikx r^ab 



(27r)^ 
1 



k^-Ml 



9^ily - (1 



a 



k^ky 



k^ - aMl 



(2.38) 
.(2.39) 



The mass M^ of the off-diagonal gluon comes from the ghost-anti-ghost condensation 
caused by the four-ghost interaction. Especially, in the SU(2) case, we obtainR 



M 



^2(^C"C") 



IGvr 



Aiie^ '''^/i^exp 



-IOtt^ 



bog^if^) 



1671 ^^^ 



(2.40) 



The SU(3) case is more complicated, see |^ for details. 



2.6 Step 5: Low-energy effective theory for diagonal fields 

We are going to calculate the VEV of the Abelian components in the given non- 
Abelian gauge theory. If the non-Abelian gauge theory can be rewritten into the 
effective Abelian gauge theory which is expressed exclusively in terms of the Abelian 



*Here we have used the minimal subtraction scheme (MS) in the dimensional regularization. 
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components only, we can calculate the above VEV in the resulting effective Abelian 
theory. 



In the previous work [O , the author has derived an effective Abelian gauge theory 
by integrating out the off-diagonal components. The resulting theory was called the 
Abelian-projected effective gauge theory (APEGT). The APEGT is expected to be 
able to describe the low-energy region of gluodynamics or QCD. In order to obtain 
APEGT, we have introduced an antisymmetric tensor field B^^ which enables us to 
derive the dual (magnetic) theory which is expected to be more efficient for describing 
the low-energy region. The magnetic theory can be obtained by the electro-magnetic 
duality transformation from the electric theory and vice versa. We have imposed the 
following duality in the tree level. 



Bl.--*{pf;. + crgr''A^^At) 



*Q 



flU 



(2.41) 



where * denotes the Hodge star operation (or duahty transformation) |3^, |4| defined 
by 



* Q,„j := -e 



fii/pa 



g^". 



First, the Yang-Mills Lagrangian is decomposed as 



£ 



YM 






4 r A^^y 



''-YM 



4 r /^^y 



(2.42) 



(2.43) 



The first piece is expanded as 



''-YM 



/^, + gf^^At 



- if 



- _f* fibcAfibAuc " j ribc Ab ac 

The simplest form satisfying the duality requirement (|2.41|) is given by 



(2.44) 



''-YM 



I-P' 



ytMu) 



2 _ }_P^gfi^J^bc^^^t>A^^^ - 



I -a' 



4 



/ {f^AlAl 



(2.45) 



On the other hand, by defining the covariant derivative with respect to the Abelian 
gauge filed, 

D,<P^ := (9^5^^ + ^/^^^A;) <?^, (2.46) 

the second piece is rewritten as 



C 



(a) 

YM 



D^Al - D^Al + gr^'AlAl 



(2.47) 



The tensor field B^i, is introduced in such a way that -B^,y-integration recovers 
the original Yang-Mills theory. Bf^i, is an auxiliary field, since it doesn't have the 
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corresponding kinetic term. However, the duality requirement leads to ambiguities 
for the identification as to what is the dual of B^iy. In fact, existence of two parameters 
p, a in ( f^.4lD refiects this ambiguity. 

In particular, when p = a, Q^^, is nothing but the diagonal component of the 
non-Abelian field strength, Q^^^ = pJ-'J^^, and hence, S^^ = ip *J^l^^- In view of this, 
the choice (|2.41|) is a generalization of that of the previous paper [^ in which two 



special cases, p = a = 1 and p = 0, a = 1 for SU{2) have been discussed as eq. (2.9) 
and eq. (2.12) respectively [0. The latter case has been first discussed by Quandt 
and Reinhardt [|T7[] . In the SU{2) case where f"''"^ = 0, the choice a = 1 completely 
eliminates the quartic self-interaction among off-diagonal gluons. 

The derivation of the APEGT was improved recently ||3^ so as to obtain the 



APEGT in the systematic way to the desired order where we have required the renor- 
malizability of the resulting effective Abelian gauge theory as a guiding principle. We 
will discuss how to choose p and a in subsection P^ . 

The strategy of deriving the APEGT is not unique. A way for obtaining the 
APEGT is to separate each field $^ into the high-energy mode ^^ and the low- 
energy mode ^^ (i.e., $"^ = $"^ + ^^) and then to integrate out the high-energy 
modes, 

~i Aa f> fii r<i na pa 

of all the fields according to the idea of the Wilsonian renormalization group (RG). 
The resulting theory will be written in terms of the low-energy modes ^^. However, we 
can neglect the low-energy modes of A" and C"', C"" due to Abelian dominance and the 
final theory can be written in terms of the low-energy modes, a* B^^, C^,C^. In other 
words, the off-diagonal components A'^ and C", C*" have only the high-energy modes. 

To one-loop level, the high-energy modes a* B^^,, C^, C^ of the diagonal components 
don't contribute to the results. Therefore, we can identify A't and C", C"" with the 
high-energy modes to be integrated out for obtaining the LEET. This strategy was 



adopted in the paper ||35|. We do not adopt this method in this paper. 



Another way is to integrate out all the massive fields 

Al^C^C. (2.48) 

Then the resultant theory will be written in terms of the massless or light fields 
a* -B* , C^jC^. The effect of the massive fields will appear only through the renor- 
malization of the resultant theory. This is an example of the decoupling theorem 
The only role of the heavy fields in the low momentum behavior of graphs with- 



out external heavy fields is their contribution to coupling constant and field-strength 
renormalization. The heavy fields effectively decouple and the low-momentum be- 
havior of the theory is described by a renormalizable Lagrangian consisting of the 
massless fields only. The decoupling theorem applies not only to theories with mass- 
less fields but in fact to any renormalizable theory with different mass scales. At 
momentum smaller compared to the larger masses, the dynamics is determined by 
the light sector of the theory. In this paper we adopt this strategy. See also section 
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(a)a^vAAAAAAA/^b,v 

P 





Figure 1: Feynman rules. (a),(b): off-diagonal propagators. The (rapidly vibrating) 
wavy line denotes the off-diagonal gluon A^ and the broken line the ghost C" or anti- 
ghost C"'. (c),(d),(e): three-point vertices, (f),(g): four-point vertices. The (slowly 
vibrating) wavy line corresponds to the diagonal gluon a* while the zig-zag line to 
the (diagonal) anti-symmetric tensor field *B^ . 



Consequently, the APEGT which was heuristically obtained in the paper 11 1 and 



improved systematically in the paper ||3^ is further modified by taking into account 
the mass of off-diagonal field components. The simplest derivation of the modified 
APEGT is to replace the massless off-diagonal propagator given in the previous pa- 
per [0 with the massive off-diagonal propagator in the QCD vacuum with ghost 
condensation.^ 

The Feynman rules are given as follows, see Fig. |T|. We enumerate only a part of 
the rules which are necessary for the renormalization at one-loop level. The two-loop 
result will be reported in a subsequent paper [R 



2.6.1 Feynman rules 

Propagators: 

^However, the following steps can be performed irrespective of the origin of the ofF-diagonal gluon 
mass. 
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(a) Off-diagonal gluon propagators: 



iD 



ab 



p'^-Ml 



9i-iu 



:i-a] 



P,^Pv 



p2 - aMX 



i;ab 



(2.49) 



(b) Off-diagonal ghost propagators:^ 



iA 



ab 



{k^ + Vd)S^'' + Voe 



,ab 



(-P - VdY + vl 

Three-point vertices: 

(c) one diagonal and two off-diagonal gluons: 



(2.50) 



i(aUp)Al{q)Al{r] 



bare 



= gr' [{q - r)^ + {r - (1 - pa)p + l}^ + {(1 - pa)p - q - '-} ^ 

where we have introduced the abbreviated notation, 

{A^ + Bp + C„\ = A^Qp^ + BpQ^p + CaQpp. 

(d) diagonal gluon, off-diagonal ghost and anti-ghost: 

i{alC\p)C\q))^^^^ = i{p + q)pgr^. 



(e) diagonal tensor and two off-diagonal gluons: 



i(*B]A''A^ 



'"' P '^/bare 



i2aglp^„^f 



iab 



where 



Iij.u,ai3 '■— -^^igfiagufs — gfi/Sgua)- 



Four-point vertices: 

(f) two diagonal gluons and two off-diagonal gluons: 



I ((/paiA^Ai)^^^^ = ig^r'^r^ 



a 



'^9^Lu9pa - 1 i9^Lp9ua + 9tia9up 



(g) two diagonal gluons, off-diagonal ghost and anti-ghost: 

-2g^nr'9,u. 



P- '^ I bare 



^.51) 
(2.52) 

(2.53) 

(2.54) 
(2.55) 



(2.56) 



(2.57) 



^ This is the ghost propagator for G — SU{2). For G = SU{5), see the paper [Q. In this paper, 
however, we don't need the exphcit form of the ghost propagator in the condensed vacuum. The 
details will be given in a forthcoming paper p7| . 
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' P+q "^N 



J,v 



r^ 








Figure 2: Vacuum polarization graphs which are necessary to obtain the APEGT. 



2.6.2 APEGT 

Thus the modified APEGT is written as[| 

^APEGT = ^ APEGT + ^{1)^APEGT, 

where C^pegt i^ ^^^ bare diagonal part, 



'-APEGT 



ifl 






7;P * KJ' 



(2.58) 



(2.59) 



and the quantum part S(^i)Capegt is obtained by calculating the vacuum polarization 
graphs in Fig. 0. In consistent with the renormalizability of the diagonal part, 



6(uC 



(I) J- APEGT 



A \J ^V) 



|*i?L.r--^(i?g2 + o 



PL] 

Ml)' 



(2.60) 



where O 



jjT ) is the modification term which is obtained in the next section. We 

use the dimensional regularization to determine 61,62 and ^3 which corresponds to 
Fig. H(a), (b), (c) respectively. Apart form the finite part, the divergent part of 61, 

^ In Fig. 0, we have omitted Feynman graphs which include external diagonal-ghost and diagonal- 
anti-ghost lines, and internal off-diagonal-ghosts (-anti-ghost) or internal o ff-dia gonal-gluon lines. 
The diagonal ghost and diagonal anti-ghost comes from the GF-I-FP term ( ^.37 ) for the diagonal 
gluon a* The higher-order terms with m ore than two external lines are suppressed by N~'^ in the 
larege N expansion, as shown in section 4T. If we neglect such higher-order terms, the bilinear 
terms in the external ghost and anti-ghost decouple fr om th e ( 2.58 ). Therefore, we don't discuss the 
contribution from diagonal ghosts and anti-ghosts to (2.5S) in this paper. 
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^2 and 6s (proportional to e ^ where e := 2 — D/2) is given by 



Si 




(2-pRaR) --H ^ — (2 - pRCTRJpRaR 

3 (47r)2 e ' 



(/i--^R)^g2(G) 
(47r)2 e ' 



52 



1 — CtR 

crR(2 - Pro-r) — o-r(1 - Pro-r) 

1+«R ^{^Ji-'g^YC2{G) 

Z (Tt 



{fi~'9KfC2{G) 
(47r)2 e '■ 



^ (4vr) 



(2.61) 



where 6^1 and 5a2 are the contributions from the graphs (al) and (a2) in Fig. || respec- 
tively, and C2{G) is a quadratic Casimir operator defined by C2{G)6^^ = f^'^D f^^D 
{C2{G) =NioTG = SU{N))!^ 

Note that the Lagrangian C^pegt is bilinear in the diagonal fields f^^, and B^^. So 
is the term S(^i)Capegt- Therefore, the divergence can be absorbed if we renormalize 
the theory. The mass generation of off-diagonal components also justifies neglecting 
the higher-derivative terms in the APEGT in the low-energy region. The additional 
term of order 0{p'^/M\) is discussed in the next step. In the SU(N) case, it is shown 
that the higher-order terms which is quartic and more in the fields can be neglected 



within the framework of the large N expansion, see section ^ 



A non-trivial consequence is that the /3-function in the APEGT has exactly the 
same form as that in the original Yang-Mills theory and is independent of the gauge 
parameter a and two parameters p, a, 



Pig) ■■= /i 



dgjf^) 
dfj, 



-h,g\^i) + 0{g' 
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-iV>0. 



(2.62) 



2.7 Step 6: Dynamical generation of the kinetic term of B 



fiv 



The antisymmetric tensor field i?^i, was introduced as an auxiliary field, since it has 
no kinetic term. In this subsection, we show that the i?^jy can have its kinetic term 
as a consequence of radiative corrections, i.e., the kinetic term of B^^ is dynamically 
generated. It turns out that the dynamical generation of the kinetic terms occurs only 
when the off- diagonal gluons are massive. This fact plays the most important role in 
deriving the dual superconductivity. 

We proceed to calculate the vacuum polarization for the B^^ field up to one-loop 
of the massive off-diagonal gluons, see Fig.0(c). Following the Feynmann rule given 
in Fig. m, the vacuum polarization of Fig. |^(c) is written in momentum space as 



n^j^,a/3(^) 



d'^p 



(2vr; 

PlP2 



^D,,,M5'''''n-'^^9n''hr^,nc 



xD.,.,(p + A;)5^i'=^| 



-2agf^'^''^I^^ 



:/3,P20-2j' 



(2.63) 



® For a special choice of the parameters, p — a — and a — Q, only the 5i has been calculated 
by Quandt and ReinhardtQ for SU{2). 
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where 



D^,{k) :-- 



k^-Ml 



QtJiu- (1 -a) 



P - aMl 



9t,v 



'^fl'^U \ f^flf^U 



(2.64) 
(2.65) 



fc2 - M2 y^" M^ J M^ k^- aM^ 
Hence the additional contribution coming from Fig. ^(c) to the APEGT is given by 

d^k 



S (I) ^ APEGT 



(27r) 



*BUk)UZUk)*BU-k). 



(2.66) 



By making use of the MS scheme in the dimensional regularization method, we 
have arrived the following expression after straightforward but tedious calculations: 
e:=2-L'/2 



^fiu,al3W 



2 " ^^v,aP 



167r' 

167r2 



^il + « 



5'n, 



u,u,a(3 I (^-^ ^ 




-x(l — x)— -75- 



In 



-^-x{l-x) — 



a + (1 — a)x — x{l — x] 



Ml 



In 



{a + (1 — a)x} 



Ml 

/i2 



^2- 
X(l -X)^r 



+ (7b - 1)(1 - a)(l -x)--iE + ln47r 



167r2 



-5'^ 



2 Ml 



[I-P),.,o.p[-]+0{e), 



(2.67) 



where Ca^*^' = fABpAB ^ jiab pab ^^^ ^^ -g Euler's constant -^e = 0.5772 ■ ■ -. See 
Appendix |C| for complete expression. 

In the neighborhood of k"^ = 0, i.e, in the low-energy region such that k'^/M'^ <^ 1, 
we use the low-energy (or large mass) expansion, e.g.. 



In 



Mj 



X{1 - X)^r 



In- 



M' 



A 



fX^ 



+ ln 



1-^(1-^)t7^ 



M' 



A/f2 °° 1 

lnl^-£-x"(l-x)" 



/^^ 



n=l 



n 



Ml 



(2.68) 



Note that this expansion is possible only when Ma 7^ 0. This is a reason why the 
dynamical generation of the kinetic term takes place. Thus we obtain 



_2CVV 
167r2 

1 k^ 



6'Hl 



^u,aj3 



e '— r^ + /o(«) + /i(«)t7^ + /2(«) 



M 



Ml 



^2M2,^^ ^^^"'"^ 



/io(a;) + hi{a] 



Ml 



O 



Mir 



(2.69) 
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where 



/o(a) 



/i(«) 



/2(«) 



1 + a^^..^ 



2 
1 



2'+(7s-ln47r-l)(l 



q; 



a^ 1 1 2 



'^)o -lE + ln47r 



/ dx x(l — x){l + ln[a + (1 — a)x]} 
Jo 



1 



a% a'^ 1 

77: < — In a 1 — 

1-a^l 3 9 9 



1 x^(l — x) 




;i + «) 



2 a + (1 — a)x 



(2.70) 
(2.71) 



— ma 1 — 

2 4 4 



11 131412,' 
a H — a a a ma 

24 3 3 24 2 



(2.72) 

(2.73) 
(2.74) 



The function / = /i(a) is a monotonically increasing function in a defined for a > 
and positive /i(a) > when a > 0. In particular, /i(0) = 0.0278, /i(l) = 1/6, and 
/i(ll/3) = 0.302. On the other hand, /2(a) is a monotonically increasing function 
defined for a > 0. Hence, < /2(a) < 1/30 for ao < a < +00 and —0.00833 < 
/2(a) < for < a < ao with ao = 0.12. Incidentally, /2(1) = 1/60, /2(2) = 0.0220, 
/2(ll/3) = 0.0258. 

In the following derivation, we don't need the explicit form of the functions hi{a). 
This is because it can be shown that the term proportional to (/ — P) (the so-called 
boundary term) does not contribute to the area law. Moreover, if we impose the 
gauge fixing condition d^ * B^^, = 0, then such terms give vanishing contributions. 
See Appendix and a subsequent paper |^ for more details. 

In the configuration space, therefore, we obtain |^ 



^(1)^ 



APEGT 



.,^C,a^g' 



<fx 



57r" 



*BUx)I^^- 






(dxd 



X\2 



Ml 



"^ /o(a)-/i(a)^ + /2(«)^^^^^ *Blp[x) 

(2.75) 

The first term in the RHS of S?^^CAPEGT gives a piece of the counter terms of the 
'^ APEGT (together with a finite part of ^3) and the remaining parts give new terms 
which can not be absorbed into the bare part, the non-renormalizable contribution. 
In the limit Ma —>■ cxd, the non-renormalizable terms disappear. It is very important 
to notice that the kinetic term for the auxiliary tensor field B^^ is generated due to 
radiative corrections and hence B^^, is regarded as the massive propagating field (at 
least in the low-energy region less than M^)- This is one of the main results of this 
paper. The implications of this fact to low-energy QCD will be discussed in what 
follows. The other contributions from Fig. ^(a),(b) can also be evaluated in the similar 

^For an antisymmetric tensor A^^, it is easy to see I^u,ai3A°'^ = A^y and A^^" I^i,^ai3 — Aap. 
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manner. Explicit form will be given in a subsequent paper, since we don't need them 
in this paper. 

Thus, we have obtained the improved version of the APEGT where the VEV of 
the functional f[a, h] of a, h is given by 



(2.76) 




using the action 



SAPEGTla, h]= I d'x [£e W + \pK~^'%,r' + C^ 



Ce[a] :-- 






4 



Ml 



^^[h] := -\{h^,^r - 2f,h^,.d,d'h^- + ^^hl^{d,d'fh^- + O (^^ 



(2.77) 



(2.78) 



(2.79) 



where we have introduced the rescaling factor K and the rescaled tensor field h^i, by 

(2.80) 
K.i^) ■■= K'/' * BUx)^ (2.81) 




and two quantities rj and 7 with mass dimension by 



27r2 



g'C2(r 



7^ := /2(a) 



.1 2n' 



g^C^a^ 



KMi. 



(2.82) 
(2.83) 



The action Sapegt[0', h] has U(l) invariance of ttfj, ^' a^ + dfj,9, i.e., 

SAPEGT[a + d9, h] = SAPEGT[a, h]. (2.84) 



Note that the Lagrangian of the modified APEGT is still bilinear in the fields, /, 



flU 



and B^^. 
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2.8 Step 7: Dual transformations 

We use an idea inspired by the (Abelian) electric-magnetic duality to calculate the 
Abelian field correlators in the cuniulant expansion (|2.25|) . Note that the Abelian di- 
agonal field correlation functions {f^^{x)f^^{y))YM can be calculated in the APEGT, 
since the APEGT is obtained by integrating out all off-diagonal components. Hence 
we have 

{fU^)fUy))yM = {fUx)fl{y)) APEGT. (2.85) 

Then the expectation value of the Wilson loop obeys 

{W{C))ym 



dficiO exp 



^ f dS^-'ix) f dS'''^{y){fl{x)fl{y))APEGT + -{'^ 



86) 



We use repeatedly the integration by parts to rewrite the expectation value of the 
electric quantity into that of the magnetic quantity as follows:|3 

ZAPEGT{F''{x)r{y)) APEGT 

Va^Vhf,^expli / d'^x {Ce[a] + Cd[h]) \ 
f^J^y^ ^expjz [d^xW'^'h J^ 



... .2KV2^^ 



/ Vafj:)hf,^expli / d^x-pK ^^'^h^.^f 

S 6 r. /■ ,, „ „„ 



j Va^Vh^, exp 1^ I d^x {^pR-^'^h^.r + C,e[ 



'2i^l/2\2 



ip 



6 

X- 



Shp^{y) 6hp^{x 
The functional derivatives are performed as 
6 6 



exph f d'^xCdlh]]. (2.87) 



exp < i / d'^xCd[h] 



Shp„{y) 6h^^{x 
JJ^^[]v[^^]h,^{x)exp[^Jd'xCd[h] 

-V[d45ix - y){5^p5,„ - 5,p5^„) + {-fV[d^]h^,{x)V[dy\hp„{y) 
xexpji /"rf^x/:d[/i]|, (2.88) 



^° Here we have used a property of the measure P/i, J'Dhjj^{- 
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where we have defined an operator, 



92 a^ 



^2 ^4 



(2.89) 



Here ^[i9] denotes quantum corrections for the h field. 

Therefore, we obtain the relationship refiecting the duality between electric and 
magnetic sector, 

{f^lu{x)fp^{y))APEGT 



+ 



'2K^nV fl-^ 



{V[dx]h,,{x)V[dx]hp,{y)) 



APEGT- 



(2.90) 



Thus the VEV of the Wilson loop in the Yang-Mills theory is rewritten in terms of 
the correlation functions of the magnetic quantity in the APEGT as 



{W{C))ym 



dficiO^W 



2J^g^p-^K f dS^'^ix) f dSP^iy) 

J Sn J Sn 



(2.91) 



' Sc •' Sc 

X {V[dx]h%{x)V[dy]h%{y)) APEGT + ■ 

This is identified as the cumulant expansion of 

{W{C))ym = et-l [dficiolexp [lip-^K^'^Jg f dS^''{x)V[d,]hl{x)\) 

•I \ L JSc J / APEGT 

(2.92) 
where [■ ■ ■] denotes the field-independent-constant part (a phase factor) in the first 
term in the RHS of ( ^.901 ) in which we have no interest. It should be remarked that 
the above result is totally independent from the explicit form of Ce[a\, even if we 
include higher order corrections of low-energy expansions. 

The above result implies that we can calculate the Wilson loop by making use of 
the magnetic theory which is written in terms of tensor field h^^ alone. By performing 
the integration over a^ field, we can obtain such a dual magnetic theory. 

By including the gauge fixing term of a^ into the action Sapegt[0', ^], 

^gfH:=-^(9X^^ 
the action SapegtV^^ h] reads 



''^J.l 1 



(2.93) 



SAPEGT[a.h] = j d^x[Cf[a]- pK-^'\'''d^h'^, + Cd[h]] , (2.94) 



where0 



Cf[a] := C,[a]+CGF[a] 



\a'p [(1 - P')9'''d' - (1 - / - p-')d^d^] < + o[j^^ 



(2.95) 



^^ It should be understood that £*°*[a] is obtained after integrating out B^C^jC^. 
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For p 7^ 1, integrating out the a^ field yields 

CM' = ^^[h] + 2[l^-'p2f K,^. {a^'d' - [1 - (1 - p')/3Wd^} d%^ 

Thus the VEV of the Wilson operator is calculated from 



{W{C))ym = Jdfic{0(exp\^2tp-'K'/'jgJ^JS^''{x)V[d,]hl 



(2.97) 



Zd' I T^Ku ^^P {' I d'^xCM'^ J[hl (2.98) 



(2.99) 



where 

J[h] := exp \2ip-^K^/^Jg f dS^"'{x)V[dMU{x) 
L J Sc 

When p = 1, Ce[a] contains only the higher-order terms. In the limit p ^ 1, the 
second term in ( p.96| ) yields the constraint 5{d^h\^) in the measure, 



{W{C))ym = Z,^ j VhlJ{d^h!^^)expi^i j d^xCd[h]]j[h]- (2-100) 

If we integrate out tensor field h^^ in SAPEGT[0',h], we will obtain the electric 
theory which is written in terms of the diagonal gluon field a^ alone. This theory can 



lead to the area law too, as will be discussed in a subsequent paper [72 . 



2.9 Step 8: Recovery of hypergauge symmetry and gauge 
fixing 

The action Sapegt[0', h] has the U(l) gauge invariance for the gauge transformation 
of the diagonal gluon field: a^ ^>- a^ + 0^6, i.e., Sapegt[0' + dO, h] = Sapegt[0', h]. In 
this section we consider the symmetry for the tensor field h. 

Now we introduce the field strength H of an antisymmetric tensor field h (the 
so-called Kalb-Ramond field), H := dh, i.e, 

H^,ux ■= dxh^^ + d^h^x + d^hx^, (2.101) 

Note that the field strength H is invariant, H^ux -^ H^^x, under the hypergauge 
transformation, h ^ h + d(, i.e., 

h>.u ^ hj,, := V + ^mC. - 5Xm- (2.102) 

We require the invariance of the measure T>h^i, under the hypergauge transformation. 
However, the Lagrangian CAPEGT[0',h] or Cd[h]' does not have the invariance under 
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the hypergauge transformation due to the existence of the mass term {h^vY- Nev- 
ertheless, we can recover the invariance by introducing new degrees of freedom^ A^ 

(2.103) 



where A^ transforms as 

A^ -. Aj := A^ 


-C 


In fact, the combination h^ := h + dA, i.e.. 




h^,. = h^u + d^K - 


-5,A 



^n 



(2.104) 



is invariant under the combined transformations, ( |2.102| ) and ( ^.103| ). Therefore, the 
Lagrangian, 

Cm[KK] := C^ih"^] = -^{h'l, + d,Al - d^KlY + y^(^;.a)' + ^(9^i^l,J^ 

(2.105) 
is also invariant, i.e., £„[/;., A] = £m[/i^, A^]. 

As we have recovered the hypergauge invariance, we need to fix the hypergauge 
invariance in quantizing the dual magnetic theory £m[^, A]. From this viewpoint, we 
adopt the condition]^ 

d'K. = 0, (2.106) 

as a gauge fixing condition for the antisymmetric tensor field.0 Under this condition, 
the derivative terms of h^^, in Cm[h^] recover the corresponding terms of Cd[h], 

{H^^xf = -^K.d^Ku - Qd,h^,.d,h'^'' - -^h\^d,d'h^'^\ (2.107) 

{d^Hi^^f = {d'dxh,, + d^d^Kx + d,d^h^,f -> /i;,(9a9^)2/i'^-. (2.108) 

Here it turns out that all the terms proportional to (1 — P) in 5%-sCapegt vanish 
under this condition. Therefore, £m[^,A] under the condition reduces to 

c^iKK] = -i(/^;^ + 9,Ai.-9.A;)2-^/^;,9,9^/^'^- + ^/.;,(9,9^)2/^'^-. (2.109) 

It is shown that the theory with Cm[h, A] reduces to the original theory given by Cd[h] 
by integrating out A field after fixing the gauge freedom of A^, see Appendix 0. 

Thus we obtain an alternative dual description of low-energy Gluodynamics in 
terms of hf^i, and A^. Especially, for G = SU{2), we obtain 

{W{C))yM = Z]^.l Jvh,,6{d''h,,) JvA,exp{tSM[h'';C]}, (2.110) 



^^This field plays the similar role to the Stiickelberg field in the massive vector theory which 
recovers the gauge invariance of the vector field. 

-"^•^In section |3.3| we discuss the relationship between this condition and the setting up of the 
previous paper |il| . 

^^In the manifestly covariant quantization of the gauge theory, we need to introduce the ghost as 
is well known. However, it is not enough for the antisymmetric tensor gauge theory, since we need to 
introduce the ghost for ghost in order to completely fix the gauge degrees of freedom. Such a theory 
is called a reducible theory. In this subsection we treat the theory in a naive manner. However, the 
result is unchanged if we take into account the reducibility of the theory. See Appendix ^. 
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where 



SM[h^; C]= f d^xC^[h^] + 2p~^K^I^Jg f dS^"'{x)V[d^]h^^{x). (2.111) 

Apart from the third term in Cm[h], the above action ( p.l05| ) coincides with the action 
of confining string proposed by Polyakov pH] in the weak field limit, see section Rl. 



2.10 Step 9: Change of variables (path-integral duality trans- 
formation) 



We proceed to rewrite the LEET ( p.lll|) into another form which is useful to derive 
the dual Ginzburg-Landau theory. 

First of all, we rewrite the surface integral in SAPEGT[h^',C] into the volume 
integral as follows. Let a = (cr^,cr^) be two-dimensional coordinate on the surface 
which is bounded by the Wilson loop C. Then 

= ^y'c?VX^'^(or)V(x(a)) 

= ]-[ rfVX^'^(a) / d^xh^,ix)5\x - xia)) 



d''xh^,{x)Q^,{x) (2.112) 

d^xh^,{x)e^,{x), (2.113) 

where we have introduced the Jacobian, 

and an antisymmetric tensor of rank two, 

e^,(x) -.= ^1 d^aXf'''{a)S\x - x{a)) = -Q.^ix). (2.115) 

We call Qfj,u{x) the vorticity tensor current which has the support on the surface 
spanned by the Wilson loop C. Note that 

Qt^-(x) = - f d^S^^^ixiam^x - x(a)), d^S''" := da^da^^j^^^^. (2.116) 
2Js o{a^,a^) 

Hence we can start from the action, 

SMlh''; Q]=J d'x {£^[/i^] + 2Jgp-'K'/'h,,{x)e,,{x)} , (2.117) 

and the expectation value of the Wilson loop is given by 

{W{C))ym = ZM[/i^;e]/ZM[/i^;0], (2.118) 
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where 



= / Vh^Jid^h^jy) exp <i d'^x 

X exp I i / d'^x2Jgp^^K^^'^h^y{x)Q^u{x) 
xJvCexp^iJd'x^ih^^^fY 

Note that the path integral transformation holds, 

-1 



1 1 



4^4 



(2.119) 



j VCf^expl^i j d'^x—ihj^^] 

J Vi^Jie^^'^^'^dpii,, - V)) exp 1^ J d^x^ii, 



where the constraint, 
is solved by 



^/liypcr 



dpii^^u - h 



liv ) 



5/^5 i.e, t^y '''flU- 



f - h = d C - d C 
Moreover, we introduce the auxiliary (Abelian) vector field 6^ by 



(2.120) 

(2.121) 
(2.122) 



die^'P^'dpii^y - h^,)) = Jvb^exp ^-tjd^x * b^,{e^, - V)| , (2.123) 

where b^^u is the (dual) field strength defined by 

b,,u ■■= d^b, - d,b^. (2.124) 

By using the identity ( |2.123| ), the integration over £^j, in ( p^.l20| ) can be performed as 
(2.120) 



Vbf, exp \i d^x* b^^{h^u) \ / 'D^fiu exp \i l d^x 



-(^p 



*b„ 



Vbp exp <{ z / d^x \-{bpyf + *b^^{h^^) 



(2.125) 



Another way of deriving the equality just derived, i.e. 



-1 



PC,exp^/rfW(/.Jj 



Vb^expi^i j d^x[-{b^,f + *b^,{h^,)\p.l2Q) 

is as follows. The argument of the exponential in the LHS is 

{h + dC,h + dC) = {h, h) + (/i, dC) + (rfC, h) + (dC, dQ ~ (/i, h) + (dC, rfC), (2.127) 
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under the condition 6h = 0. The last term decouples after the Gaussian integration 
of (. On the other hand, the argument of the the exponential in the RHS is cast into 



d X 



ib 



/.w ) 



* V( V) = (db, dh) - {*db, h) = (b, Sdb) - (6, *dh). (2.128) 



Suppose the Lorentz type gauge condition 6b = 0. We introduce the NL (zero- form) 
field 0. Then the Gaussian integration over b^ field yields 



(6, Ab) - (6, *dh) - {6b, 0) = (6, Ab) - {b, *dh + d(t)) 

1 ,, , ,\ /, ^d^ , , . 6d 
*dh + d(j)) = {h,—h) + 



(*dh + 



A 



A 



A 



0)~(/i,/i) + (0,0),(2.129) 



under the condition 6h = 0. In this derivation, we must insert the constraint 6{d^b^) 
in the measure "Db^. The identity implies that there are many ways of extracting the 
transverse modes of h. 

Thus, the theory is rewritten in terms of b^ and h^y as 



vh^Jid^h^y) 



d'^x 



*bf,uhf,u + -z^^Hx^uf + —l{^^Hx^,yf 



ZM[b,h-Q] 

Pfe^exp*}? / d'^x \-{b^y) 

X exp < i 

xexp jz f d^x2Jgp-^K^/^h^y{x)e^y{x)\ . (2.130) 

By change of variable b^„ -^ b^^ — 2p^^K^^'^Jg * Q^u, we arrive at the expression, 
ZM[b,h;Q] 

Vb^exp \t fd^x \-{b^, - 2p-^K^'^Jg * e^,)^^ 



X / Vh^Jid^h^y) expli d'^x 



*b^yll^y + 



Y^iHx^uf + —^{d^Hx^^f 



(2.131) 



2.11 Step 10: Dual Ginzburg-Landau theory in the London 
Umit 

First of all, we examine a special case 7 = cxd for simplicity, although our derivation 



suggest 7 < 00. So the last term in (|2.131|) is neglected. The case of a finite 7 is 
treated in the next section. We change the variable h^^ into the new variable V^ as 
follows. ^^ 



T>h^y6{d''hf,y) exp < i / d'^x 



12r]- 



■{H, 



Xfii/ ) 



+ *bayh 



/lU'"/!!/ 



^5 From (5/i(^) = 0, there exists a three-form F^^) ^^^-^ ^^j^a,t /i'^' = SY^^^ =6* W'^^^ = *dW^'^\ 
Then ff(3) _ ^^(2) ^ d^dW'-^^ = *6dW^^'> = *V^^\ or l/(i) = *H^^l Therefore, 6V^^'> = S*H^^^ = 
^dH^^'' = *dd/i(2) = 0. 



28 



Vh^y5{d''hf,y) exp < i / d^x 



12r]- 



:{Hx^uf + e>''"'%,d,h,^ 



VV^5{d^V^) expli (fx 



^X + ^b,V^ 



(2.132) 
(2.133) 



since the constraint dfj_V^ = can be solved by an antisymmetric tensor field in the 
form, 



yt^ := -e'"'f"'d^h 



iivpcj 



H,. 



d, * h^r 



(2.134) 



^ — (X ■ " jju ^ — — '^ U pU 

The massive antisymmetric tensor field h^y denotes the massive spin-1 field V^ whose 
canonical mass dimensions is three.[3 In this step, the number of independent degrees 
of freedom is conserved, since V^ and h^y have three independent components. 

Furthermore, the integration over V^ is performed after introducing the new vari- 



able 9 to remove the delta function of the constraint d^V^ 



0,13 



^2.133) 



VV^, VOexpl / d^x 



VVf, I Veexp\i (Tx 
1 






27^2 ^ 



r]\2b^ - d^ey 



= fvOexph fd'^x 
Finally, we obtain the dual Abelian gauge theory, 

ZM[b,e;e] 

= J Vb^ I V9 exp h I d^x --{b^y + t 
where 



(2.135) 
(2.136) 
(2.137) 



S N2 



lv'ib,-d,i 



d'^hl^x) 



-^K^'^vmJ: 



p ^K^'-v[a\j;, J^ 



btix) := Ap-^K^^^Jg*Q 
■1 dxuir 



' iiu 



flU 



X) 



:2.138) 



(2.139) 



Ug I dr^^S^ix - a; (r)), (2. 140) 
dr 



and g is the Yang-Mills coupling constant of the original Yang- Mills theory. 
This model has dual U(l) symmetry, say f/(l)m symmetry. 



(2.141) 



This model is identified with the London limit X —>■ oo oi the dual Abelian Higgs 
model or the dual Ginzburg-Landau theory with the Lagrangian, 



-1 



2\2 



CDGL[b, 0] = —ib,y + b%Y + \{d,- tgmbM - A(|0r - v') 



(2.142) 



^^The massless antisymmetric tensor field stands for the massless spin-0 field, see 

^^ In view of the definition ( ^.134 ), the V^ plays the similar role to the magnetic monopole current 



k^ defined in the next section. Here the variable 9 is introduced to keep the constraint 9'^fc^ = 0. In 
other words, the introduction of keeps the magnetic C/(l)m symmetry. So, putting 9 — breaks 
the U{l)m symmetry. 
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where Qm is the magnetic charge subject to the Dirac quantization condition, 

9m9 = 47r. (2.143) 

The London hmit is equivalent to putting \4>{x)\ = v = const., i.e, (l){x) = vexp[i9{x)]. 
The dual U(l) symmetry is broken in the London limit, 



(2.144) 



This corresponds to the infinitesimaly thin flux tube connecting the quark and anti- 
quark. In the London limit, the Higgs mass m^ = 2^/Xv diverges, i.e., m^ = oo or 
mT^ = 0. This is the extreme case of the type II superconductor where m,^ > rrib. It 
turns out that the mass m^ of the dual gauge field b^ is given by rj, 




/K V2471 

9 



(2.145) 



The monopole condensation is shown to occur in section ^ 



3 Final step: Dual Ginzburg-Landau theory of the 
general type 



In this section, we discuss how the LEET given by ( p.l31|) is related to the dual 
Ginzburg-Landau theory of type II. In section |^, we give another evidence of equiva- 
lence between the LEET ( |2.131| ) and the dual Ginzburg-Landau theory on the border 
between type II and type I. 



3.1 Dual gauge theory 

We return to eq.( p.l3lD : 
ZM[b,h;Q] 



X / Vh^.Jid^h^^) expli d^x 



*b,.h,^ + j^,iH,,^r + ^,id'H, ^^ 



l\fj,u) 



By making use of the change of variable ( 2.134 ), we change the variable hf^^y into the 



new variable V, 



tj.1 



/ Vh^y5{d''h^y) exp < i / d^x 
/ Vh^y5{d''h^y) exp < i / d'^x 
jvV^5{d^V^)expUjd*x 






2b^V^ 



1 -.2 1 



^x-^.^^y^-d.v,? 



(3.1) 
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After introducing the Lagrange multiplier field 9 for the constraint dyV^ = 0, the V^ 
integration is performed as 



(3.1) 

fvOexpli f d'^x 



^ "2 1 /Q TA i^ T/ \2 



(26, - 0,6)^^ - —V; - ^(9,K - d.V,) 



l^'^^-'^'^A^^^^ - 



gi,^ _ !l_g^g^ ] (26^ - dj 



(3.2) 



Then we obtain 



ZM[b,e;e] 

Vb^Sid^b^) fveexpli f d^x 



7(V + &; 



S N2 



7 



+ 2^^'^ - ^m^) a_;4/^2 1^^^'^ - fi^^^^ 1 (b. - 5.^ 



?7 

7 



(3.3) 



where we have inserted the delta function 6{d^b^) for fixing the gauge for 6,. Note 
that ( |3.3| ) reproduces the London limit when 7 -^ 00. In the case of finite 7, 



l(6,-a,%Ml-^A 



-1 



7 



^^"-^9^5'^j(6.-9,f 



2 ^ ' I 74 



1 



6^ - -r^^^A^ + 0(A2) 



§^'&/.^^ + ^^^^A^^ - Iv'OAe + 0(A2) 



(3.4) 
(3.5) 
(3.6) 



where we have used d'^b^ = 0. Thus we arrive at a LEET of the Yang-Mills theory, 



^APEGT 



[b,e;Q] 



Vb^5{d%^) expli d^xCK[b, 6] 



^K[bA := -^ (1 + ^Vm' + ^ A&" - ^^'^A^- 



(3.7) 
(3.8) 



The LEET just obtained is of the same form as ( |2.138| ), except for the renormalization 
of the kinetic term of the dual gauge field. The dual gauge field becomes massive, 
whereas the 9 field remains massless. This is reasonable, since the field 9 corresponds 
to the Nambu-Goldstone (NG) mode associated with the spontaneous breakdown of 
the magnetic U(l) symmetry. 



3.2 Low-energy effective theory of dual Abelian Higgs model 



In the following we discuss how this theory (|3.8| ) is related to the dual Ginzburg- 
Landau theory of type II. We remember that the London limit itlh -^ 00 corresponds 
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to the limit 7 = oo. Therefore, we expect that the LEET with the Lagrangian ( p.8| ) 
can be reproduced from the dual Abelian Higgs model, in the region = m^ -C rrih -C 
rriH- 

We show that the dual Abelian Higgs model with the Lagrangian, 

CDGLlb, 0] = ^{b,. + b%f + \{d, - igmbM" - A(I0I' - v^)\ (3.9) 



reduces to the LEET with the Lagrangian (|3.8| ) in the low-energy region. We adopt 
the renormalizable gauge PTl E3], 



5% + e"^6¥'2 = 0, (3.10) 

where the scalar field is parameterized as 

(t>{x) = -^[v + ^i{x) + i^2{x)\. (3.11) 

The GF+FP term of the renormalizable gauge is given by 

C-GF+FP = -iz^{d^b^ + ^mbip2f + ic{d^ + ^ml)c + igm^nibCCipi. (3.12) 

Even in th Abelian gauge theory, the renormalizable gauge require the FP ghost and 
anti-ghost which have a non-trivial interaction with the Higgs scalar ipi. A merit of 
the renormalizable gauge is that the mixing term inhb^d^ip2 between b^ and (^2 in 
Cgf+fp cancels the same term in the original Lagrangian, 

\DAbW = \{d>.-t9mb^W (3.13) 

= 2*^^'^'^! + 9mbf,(p2y + -^{d^ip2 - 9mbf,ipif 

2 2 

-gmvb^id^Vi + 9b^.^l) + ^^m^'^- (3-14) 

Note that ip2 is the would-be Nambu-Goldstone (NG) boson which corresponds to 
the phase factor 9 in the polar coordinate, 

0(x) = i=[t; + p(x)]e^^(-). (3.15) 

However, we don't use this parameterization in this section, since it is not suitable 
for performing the loop calculation. Then the total Lagrangian of the DGL theory 

CDGL[b, 0, c, c] = ^{b^,f + I (9^ - igmbM"" - ^(l^r - ^')' + 'I^gf+fp, (3.16) 
is decomposed into the free part and the interaction part as 

-\{d^b, - d,b,)' + \mlb,b^ - ^(a^&,J^ + i-c{d' + iml)c, (3.17) 
C-i = gmbf,{d>'^i^2 - '9^<^2<^i) + g^vb^b''^, + -glb^b^'iA + A) 

-Xvifiiiipj + ifil) - -{ipl + iflf + igm^rribcapi, (3.18) 
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where the dual gauge boson mass mf,, the Higgs scalar mass m^^ and the would-be 
Nambu-Goldstone mass m^^ are given by 

rrih := gmV, rri^^^ = 2/i^ = Af^, m^^ '■= C'^l- (3.19) 

The Feynman rules is given in Fig. ^. The propagators are given as follows. 
Propagators: 

(a) Higgs scalar ipi propagator: 

iDi{k) = TT ^ -, (m 



j^2 _ ^2 ^ _|_ ^g ' 



^1 



2/^2 = Xv^). 



(3.20) 



(b) Would-be Nambu-Goldstone boson ip2 propagator: 



iD^ik) 



f^2 _ ^2^ _^ ^g 



rn^2 =^^b =^9m" )■ 



(3.21) 



(c) Gauge boson h^ propagator: 



iDf^uik) 



k'^ — ml + ie 



9txv 







kfiky 



k^ — ^ml + ie 



(3.22) 



(d) Ghost (anti-ghost) c (c) propagator: 



iD,{k) 



k'^ — C.ml + ie 



(3.23) 



The relevant vertices in the Feynmann rules are given in (e) to (j) of Fig.|^. 

The basic strategy of showing the equivalence is to integrate out the massive 
scalar (dual Higgs) field ipi with mass squared m^^ = 2/i^ = \v^, since the theory 



is written in terms of b^ and 9 



^2/v)- The renormalized mass rriH 



m, 



1^2 



of the heavy field is made large while all other parameters are held finite. The irif, 



and m^^ are the masses of the light fields. The decoupling theorem [18| asserts that 
phenomena on energy scales much less than the Higgs mass mn = rn^p^ are described 
by a low-energy effective theory with the Lagrangian 



u 



-\{d^K-d,h,f + KnXh^ 



Z,{dX - dub^r + -mlZ,hy^, 



(3.24) 



where we have substituted the renormalization relation, h,, := Zj b^. On the other 



) "M 



hand, the renormalized Lagrangian with the counter term is given by 



C, = --idX-d.h 



,. •krn^)Xb'''-W^b"-9^b".^" + Wb^.b'''- (3.25) 
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Figure 3: Feynman rules for the dual Ginzburg-Landau model: Propagators: (a) 
Higgs scalar ipi propagator, (b) Would-be Nambu-Goldstone boson ip2 propagator: 
(c) Gauge boson 6^ propagator, (d) ghost C propagator; Vertices: (e) ipi — ip2 — b, (f ) 
ipi-b-b, (g) (fi - ipi - b - b, ip2 - f2 - b - b, (h) ipi - ipi - (fi, ipi - ip2 - <^2, (i) 
ifi-(fi-(fi~ If I, y:>i-(fi-(f2- f2, (f2-f2-f2- f2, ii) fi - c - c . 



Equating ( p. 241 ) and ( |3.25| ), we obtain the relationship. 



Sh := Zh 



5r, 



Zbtnl 



[771 
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We use the renormalized perturbation theory in the (dual) coupling constant 

47r 

9m ■= ■ 



(3.26) 



(3.27) 



In order to see coincidence of the resultant theory with Lagrangian (p.8|), it is 
useful to choose the Landau gauge ^ = in which the would-be NG boson ip2 is 
massless, since m^^ = ^m^. In Appendix ^ we evaluate one-loop graphs for the 
self-energy of the would-be NG boson ip2 and the vacuum polarization of the dual 
gauge boson 6^ as shown in Fig. ^ and Fig.^ respectively, based on the Feynman 
rules given in Fig. ^. The radiative correction causes the mass renormalization and 
wavefunction renormalization of b^ field, while the would-be NG boson field ip2 field 
remains massless for ^ = 0. 

Now the renormalized mass is written as 



[771 



R\2 



Zb7nl - 6m = ml + ml6b - S„ 



(3.28) 
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where {jnfY is of order g'^ml, whereas 6m is of order g^^ or g'^m'^j. Actual calcu- 
lations in Appendix ^ give 



Hva 



N, 13 



(47r)2 \ 3 



+ — + 6 In 
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m 



fi- 



+ 0{gt), 



%: 



(47r 



™m2 



;^'"b 



'iN, + ^] 



171% In ^ — m? In ■ 



m 



H 



mt 



8 m? 



21n- 
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H 



/i^ 



(47r; 



2"^/, 



where 



-N, + - + In 



N. 



m 



H 



1 



+ 0{g: 



m/i 



+ ln47r — 7£;. 



(3.29) 
(3.30) 



Substituting (|3.29|) into ( |3.28|) , we obtain 
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(47r; 



2m^ 



A^e + ^777 + 5 In 



m 



H 



36 



/i- 



In 



mt 



/^^ 



9m 2 

74^0^ ""' 



-(A^. + l)- 



m^ In ^^ — ml In ■ 

2 2 

m^jj — mil 



+ 0{g'J. (3.31) 



In order to keep the physical mass mf of 6^ finite even for the large dual Higgs mass, 
m,H ^ 00, we must let m^ have a term proportional to g^^m^j^, 



m. 



finite term + 



9™ 9/4,, 49 ^, m?, , m? 



-m,i 



(4vrp"^^ ^3 
Thus we arrive at the conclusion, 

(m^ )^ = finite term in m,b 
4q2 

i/rr 



^■ + s + ^'"7? 



-In- 



/i^ 



(3.32) 



:m. 



9 1 "itr 9 1 "*f 

3_ _ m|^ In -f - m^ In -4 



^(iVe + l) 



(47r)2"^'' 
The comparison with ( p. 81 ) implies that 

/I 



m|^ — m^ 



+ 0(^^^). (3.33) 



1 + ^, (m,)^Zfe = r]'. 



r 



3.3 Comparison with the previous work 



(3.34) 



In the previous paper ||Tl[, we have used the Hodge decomposition for the two- form 
B, 

B = db + *dx, (3.35) 

with two one-forms b and x- Then the definition of h yields 



h := *B = *db + **dx = 5*b — dx-, 



(3.36) 



35 



which leads to 

Sh = -Sdx, (3.37) 

H ■.= dh = d6*b= *5db, (3.38) 

where we have used dd = = 66. Moreover, 

6H = 6dh = 6* 6db = *d6db = *d{6d + d6)b = *dAb. (3.39) 

Hence, if we require 6h = 0, then A^ = under the gauge fixing condition 6x = 0. It 
is known |3^ that a p-forni uj is harmonic {Auj = 0) if and only if uj is closed [du = 0) 
and co-closed {600 = 0). The harmonic form u does not exist on the topologically 
trivial manifold, since the dimension of the set of exact p-forms is equal to the Betti 
number, i.e., dim Harm^(M) = b^. In this case, x is divergenceless and rotation free 
vector field on the four-dimensional Minkowski spacetime. Hence we can eliminate 
the variable x cind hence x does not appear in the result. This corresponds to the 



situation discussed in the paper |11 



It is instructive to compare the above result with the previous one. Replacing 
V^ = du * h^^ with the magnetic monopole current A;^ (this definition is suggested 
from h ~ *B since k ~ 6B), we obtain the theory. 



ZAPEGT[b, fc; 6] = / Vbf, I Vkf, expl i J d'^x 

+b,k>^ + ^kl + ^{k,Ak) }. (3.40) 






After integrating out the dual gauge field 6^, this action leads to the theory of mag- 
netic monopole written in terms only of the monopole current fc^. See section ^. 

3.4 How to determine the parameters p, a, a 

It has been shown [|l7i [TT| , |35| that the renormalization of the Yang-Mills coupling 
constant g does not depend on p, a and the gauge parameter a. Therefore, the /3- 
function is also independent of the choice of p, a and of the gauge parameter a. The 
resulting /3-function exactly coincides with the one-loop /3-function of the original 
SU{N) Yang-Mills theory 

(3{9K):=fi^ = -9KP^lnZ, = -j^^gi + 0{gl), bo = ^C,{G), (3.41) 

exhibiting the asymptotic freedom. Moreover, some of the anomalous dimensions have 
been calculated.f^ The anomalous dimensions of the fields in the SU{N) Yang-Mills 

^^The anomalous dimensions of the fields a* i?* and the parameters p, /3 are calculated in the 
previous paper pq]. To obtain the anomalous dimension for a, a, we need to calculate more Feynman 
graphs than those in [B5|, see p3|. For G — SU{2), Za and Zc; were calculated by Hata and Niigata 
|3^ and Schaden |19 , Za by Schaden [Q and Zc^ by Quandt and Reinhardt [|l7|] . 
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theory are evaluated as 

la^ig) ■- 

lA-ig) ■- 
ic-ig) ■- 
ic-{g) ■■ 



^.|ln^.^^C.(G)(|^ [sum. 






22 _ 9 an 

y ^ 2 ~ T 



[SU{N)], 



-/3fi)c2(G)-^ [SU{2)]. 



(Att) 



(An) 



[SUi2)], 



2W '' 



(3.42) 



where ? denotes that the result is not yet unavailable. For the parameters p,a,a,P 
and C, 



ipig) 



dpR ^ ^ 7 

ii _ ^ + 2— - ^~^^ (^^ - '^ 
6 2 PR 2 ' 



7a (^) 

7a (^) 
7/3 (^) 

7c (^) 



-pR 

dan 



.PR 



C2(G') 



[All) 



[SU{N)l 



dp 

dag 

dfj, 

OPr 

djj, 



-crRfi--\nZ„ =?, 

-c^rPtt- In Za = -an \-6- — + aR] —^ 

OjJ, \OiR 6 J STT^ 



[SU{2)1 



-27a«(^)/9R [SU{N)l 



-PrPtt- In Z^^ 

OfJ, 

d / 3 



dfi 



a 



[SUi2)]. 



(3.43) 



-R/ 



Note that the parameters p, a, a, /5 do run and changes according to the scale p.. We 
would like to obtain the renormalization scale p independent result. The simplest 
way is to search for the fixed point for these parameters on which the parameters are 
kept fixed irrespective of the renormalization scale p. The fixed point is determined 
by solving the simultaneous equations, 7p(5') = 0,7o-(5') = 0,7^(5') = and 7^(5') = 0. 
Finally, we give a simple argument how to determine the parameters p, o", a in our 
theory. The parameter p obeys the differential equation. 



ipig) ■= p 



dpR 



11 1 + a o\ 

y + -^—^ j Pi? - (3 + a)a 



C2iG) 



(47r) 



This implies that the point p = p^j 

p = p 



a 



3 "•" 2 "^ 



(3.44) 



(3.45) 
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Figure 4: Higher order terms for obtaining corrections of APEGT. 



is the infrared fixed point, as far as y + ^-^-a"^ > 0. At least in the present stage of 
investigations, the simplest choice of the parameters is 



aR 



c^R = 1- 



(3.46) 



The choice a = 1 greatly simplifies the evaluation and the expression of the result for 
the vacuum polarization of 3^,^. The choice a = 1 completely eliminates the quartic 
gluon interaction, see ( |2.45|) . Usually, the auxiliary field B^j^i, is introduced so as to 
achieve this situation. In this case, the value 



6 

PR= ^ 



(3.47) 



is the infrared fixed point. However, ( p.43|) implies that aR^n) monotonically decreas- 
ing in /i and that there is no fixed point for a in the SU{2) case, although the SU{N) 
case can be different from the SU{2) case. 

The complete list of the anomalous dimensions in the SU{N) case and the details 



of the RG properties of the APEGT will be given in a subsequent paper |72 . 



4 Estimation of neglected higher-order terms 

We proceed to discuss the issue whether the neglected terms in the above derivation 
do not invalidate the above result in the range of parameters and the energy scale in 
question. 
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4.1 Higher-order cumulants and large N suppression 

In this paper we have obtained the APEGT which is bihnear with respect to the 
diagonal fields, a^ and B^^. Hence the bilocal approximation in section |2.3| is exact 
within this framework. In the framework of large A^ expansion, moreover, neglecting 
higher-order cumulants can be justified as follows.[^ In the large N expansion, A : = 
g'^N is kept fixed. 

For example, we consider the diagram in Fig.^ with n external B^y lines. Accord- 
ing to the Feynman rules in Fig. |1], it corresponds to 



nJ (27r) 

')diC2 ( _ 

cnp2\t^ ^ "'^fy ^y^ J -^ t^2^2,P20-2 



xDi\';l{p + h)[-2gaP''''H, 



xDt'TpliP + ki + --- + fcn-i)[-2(7af-"'^" /,„,„,,„.„]. (4.1) 

This quantity is proportional to the factor, 

^nSjdnCl rilCldl SjdlC2 £i2C2d2 . . . rd„_lCn £inC„dn (A 2] 

For n = 2, the factor reads 

g2^d2Cipicidi^diC2p2C2d2 ^ g2 piciC2 p2C2Ci ^ _g2(J^^ili2 ^ _g'i]^^hi2 ^ 0{N°). (4.3) 

For n = 3, it is identically zero, 

3 rdsci riicidi ^diC2 p2C2d2 rd2C3 pac^ds __ 3 piciC2 H2C2C3 p^csci _ q f4 4) 

For n = 4, it is of order 0(A^~^) for large A^, since 

4 rd4Ci picidi ^diC2 p2C2d2 xd2Cs pac^d^ ^d^c^ p^Cidi 
__ 4 piCiC2 p2C2C3 p:jC3C4, pACAC\ 

= C;4^^(^^ni2^i3i4 ^ §ili3§i2i4 ^ §ili4§i2i3^ _ ©(A^'^). (4.5) 

The last equality is obtained as follows. From the symmetry under the exchange of 
the indices ii, i2, "^3, "^4, we can put 

piC\C2 P2C2C3 P3C3C4 «4C4Ci _ /I /'rnJ2X«3*4 _|_ Xn«3X«2M _|_ XnMX*2«3'\ /A Q\ 

By contracting 23 and 24, we obtain 

r tJC _ piC\C2 P2C2C3 P3C3C4 P3C4C1 __ _ piCiC2 P2C2C3 rciC3 _ T\T^iii2 fA Y] 

RHS = A[6'''''{N-l) + 25'''''] = A{N + 1)5'''\ (4.8) 



The author would hke to thank Toru Shinohara for helpful discussion on this point. 
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where we have used the formula ( [A.3| ). Hence we obtain A = N/{N + 1). This 
argument can be extended to arbitrary n. Thus it turns out that the higher-order 
terms with n external lines of the tensor fields B are suppressed by 1/A^^ in the large 
A^ expansion for n > 4. 

Thus the contribution from the diagrams with n external lines of diagonal fields is 
suppressed by a factor l/N"^ for n > 4 where n is the total number of external diagonal 
gluon fields a^ and external tensor gauge fields i?^^ in the off-diagonal one-gluon- 
loop or one-ghost-loop diagram. This is because every three-point vertex (c),(d),(e) 
in Feynman rules has a common factor gf^"'^. In the leading order of the large N 
expansion, we have only to take into account the diagrams with two external lines 
and hence the resulting APEGT is bilinear in the diagonal fields, a* (or /* ) and 
S* . In this limit, therefore, the bifocal approximation of neglecting the higher-order 
cumulants is exact within our approach. In this sense, the bilocal approximation is 
consistent within the framework of the APEGT. This situation is in sharp contrast to 
the bilocal approximation in the analytical treatment of the stochastic vacuum model. 



although the validity is confirmed by the numerical calculations on a lattice ||30 



4.2 Higher-order terms of low-energy or large mass expan- 
sion and the decoupling theorem 

We have neglected higher-order terms of the large mass or the derivative expansion in 
powers of j9^/M^. This approximation will be valid in the low-energy region below Ma- 



This is considered as an example of the Appelquist-Carazzone decoupling theorem 1 18 

First, we recall the case of QED. A typical example of applying the decoupling 
theorem is a derivative expansion of the photon effective action (known as the Euler- 
Heisenberg Lagrangian) obtained by integrating out the electron field in QED as 

= -i\n [dtp] [d(f)] exp {iSqed} 



^ 90(47r)2M4 ' " ^ 



3(477) 




.2 \3 



where M is the electron mass and 

2 M^ 

z(/i) := - + ln47r-7s-ln^. (4.11) 

c Lt 

It turns out that the terms that do not decouple in the M —>■ oo limit have the 
same form as those appearing in the original Lagrangian and therefore they can be 
absorbed in the wavefunction renormalization. The new structures appear as non- 
renormalizable terms and they vanish in the limit M ^ oo. This is an example of 
the decoupling theorem |jT8|. The theorem states that, under some given conditions. 



the effects of the heavy particle only appear in the light particle physics through 
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corrections proportional to a negative power of M or through renormahzation. The 
vahdity of this approach is hmited to energies much lower than the mass M. 

In the low-energy effective action of the Yang- Mills theory, all the terms including 
more than two diagonal fields can be suppressed in the large N limit, as we discussed 
in the above. Of course, this argument does not hold for the relatively small A^. Even 
in the case of not-so- large A^ (e.g., iV = 2), however, such terms are suppressed by 
the power of the inverse (off-diagonal-gluon ) mass l/M^. In the limit M^ — »■ c>o, the 
off-diagonal gluons affect the physics described by the diagonal gluons only through 
renormahzation. In the case of small A^, the validity of our approach is limited to low 
energies lower than the off-diagonal gluon mass Ma- More quantitative estimate of 
the higher order terms will be given in a subsequent paper [^ . 



5 Magnetic monopole condensation and area law 

In this section, we show that the LEET given by ( p. 1311) reproduces the same results 
as those obtained by the supposed dual Ginzburg-Landau theory of type II. Here 
the dual Ginzburg-Landau theory can include the type II on the border of type I. 
Therefore, two extreme limits, London limit and Bogomol'nyi limit, are special cases. 

5.1 Monopole action and monopole condensation 

In order to obtain the monopole action, we return to the action (|2.117|) , 



Suih; 6] = jd^x [Cm[h] + 2Jgp'^K^I^h^,{x)Q,,{x)} , (5.1) 
where the expectation value of the Wilson loop is evaluated as 

{W{C))ym = ZM[h;e]/ZM[h;0], (5.2) 

ZM[h;e] := J Vh^J{d''h^,)exp{iSM[h;e]} . (5.3) 

We define the monopole current A;^ by 

K ■= amd" * V = grnK'/^d'B,,,, (5.4) 

which satisfies the (topological) conservation law, (9^A;^ = 0. Conversely, the tensor 
field h is written in terms of the monopole current k as 

'^fii' = 9m ^fJ.upa-T^'J "^ = 9m 7: ^ l-W pa 1^ [O K —OK), (5.5) 

which is subject to the constraint, d^h^y = 0. It is easy to show that 

2 

/ d^x^Qi^yf = (k, A-^k) + (Sh, A-^h) ~ (k, A-^k), (5.6) 

-Jd'x^{H,,,r = {k,k), (5.7) 

2 

I d^^9^(d'H,^,f = ik,Ak), (5.8) 
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where we have defined A := d5 + 5d and used 5h = 0. Thus, we obtain the magnetic 
monopole theory as a LEET of Yang-Mills theory, 




where t] and 7 are given by ( 2.82 ) and ( 2.83|) respectively. Hence the expectation 
value of the Wilson loop is obtained from 



{W{C))ym = ZMp[k; E]/ZMp[k; 0], 



(5.10) 



and 




(5.11) 



where S^ denotes the four- dimensional solid angle under which the surface^ S with 
the two-dimensional coordinate (a^, a^) is seen by an observer at the point x, 



^e^upad:J^dS'''{x{a))A-\x - x{a)) 



1 



3 o '-fj.upa ^ 
57r^ JS 



d: / dS'-'ixia)) 



1 



(x — x(o"))^ 



(5.12) 
(5.13) 



Here note that the total number of independent degrees of freedom is unchanged 
under the change of variables from h^y to k^ and that the associated Jacobian in the 
integration measure "Dk^ is field independent and hence omitted. 

From the monopole action, we can demonstrate that the monopole condensation 
does really occur in the sense {k^kp) 7^ 0. In the London limit 7"^ = 0, the propagator 
of the monopole current is given by 



{kp{x)ky{y)) =g^,y 






v%. I -^.J^e^'^^-y^ 



(27r)4z 7^2 - p2 



(5.14) 

Therefore, we obtain non-vanishing monopole condensation for non-vanishing off- 
diagonal gluon mass. 



{kf,{x)kp{x)) 



d'^p rj^p^ 1 

(27r)4r/2-p2 = ^ 



7] 



(ln47r — 7^; + 1 — ln?7' 



(5.15) 



where we have used the MS scheme of the dimensional regular izat ion. This should be 
compared with the mass m^ of the dual gauge field h^. An close relationship between 
the monopole condensation and the mass of the dual gauge field was conjectured in 

^°An apparent S dependence should drop out after summing over branches of the multivalued 
potential. 
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the previous work |]rT|. In fact, the above propagator for the monopole current leads 
to 



{k^{x)K{y)) = -7f6^ 



^,.K-j'^uyyj/ ; -/.. , ^^^y 






*(-») = _^^i^j(x_y)+o{,,'), (5.16) 



which is nothing but eq.(4.22) predicted in the previous paper [|TT|. In this paper we 
have shown that the origin of monopole condensation is the existence of off-diagonal 
gluon mass Ma which provides also the mass rrib of the dual gauge field b^, 



(k^k^) ^ ^ Ma ^ ^ mfe ^ 0. 



The monopole action ( |5.9| ) is written in the form, 

1 



SMp[k] 



H 



{k,DUp')k), 



where the inverse of Dm{p) is the propagator given by 

p2 p 



D;n\p') ■■= { i, - -, + ^] ' ' ^' 



j)'^ jl'^ '■y^ 



X 1 9 9 9 9 i ' 



where 




X 



1 1 

2 2 

m| — ?Ti2 



r] 



1 - V/74 



Finally, the monopole condensation is calculated according to 
which yields 



{kuk^} 



X 
47r2 



(In 47r — 7^; + 1) (m^ — ■022) — rn^ In ml + m\ In mi 



Thus we obtain non-zero monopole condensate. 



(5.17) 



(5.1^ 



(5.19) 



(5.20) 
(5.21) 



(5.22) 



(5.23) 
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5.2 Area law of the Wilson loop 

Thanks to the NAST, the Wilson loop operator has an alternative form expressed 
in terms of the monopole current. The expectation value is expressed by ( ^.11|) . 
By making use of the monopole action ( p.l8D , the VEV of the large Wilson loop is 
calculated by performing the Gaussian integration over k^ as 

{W{C))yM = exp i^-^-i2Jgp-'K'/'nE„ Z^^^S^)} . (5.24) 

It is not difficult to show (see Appendix ^ that ( |5.24|) leads to the area law, 



{W{C))YM^exp[-astA{C)], 



(5.25) 



for the large Wilson loop. The string tension is obtained as 



0"st 



J'g' - 



m. 



2tt 



p-'Kx\n'-^, 



rrio 



(5.26) 



This is one of main results of this paper. The static potential is defined for the 
rectangular loop with side lengths T, R as 



V{R):=-\im -\n{W{C))yM. 

1 ^OO J 



Then we obtain the hnear static potential, 



V{R) = astR, 



(5.27) 



(5.28) 



for large separation R. This result is consistent with the claim that the QCD vacuum 
is the dual superconductor (of type II). This result indicates that quark in any rep- 
resentation is confined in the SU(2) case. For SU(3), our result can be applied only 
to quark in the fundamental representation due to a restriction coming from NAST. 
The factor p~^K is estimated as follows. In particular, when a = 1, /q is given by 



/o 



In- 



Ml 



7e + In 47r 



In 



9' 



1 + 



167r2 



(5.29) 



where we have used the expression of off-diagonal gluon mass ( p.40| ) . Then K reads 



K -.= 1 



27r2 



■^Vo 



24 

1 + — a 
11 



' ^^ a'(ln.^ + l\:=K{g), 



27r2 



levr^ 



(5.30) 



since C2 = N and 60 = ^^ for G = SU{N). For the choice a = 1 and p = 6/7 in 



( |3.47| ), the factor reads 

p-^K 



49 
36 



35 Ng^ / g^ 
11 27r2 I "^ 167r2 



(5.31) 
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which is positive for < g < 9.6 and monotonically increases from [4.3, 6.5] for 
g e [0,4.6] with a peak aX g = 4.6 and monotonically decreases with [6.5,0] for 
g G [4.6,9.6]. Therefore, we can put p~^K = 5 for gf = 2 in the SU(2) case. The 
numerical estimation of x, ^^i, ^^2 will be given in section ^. 

From the viewpoint of dual Ginzburg-Landau theory, two constants mi and m2 
may be regarded as the coherence length m^ and the penetration depth rrih (apart 
from a factor \/2). Therefore, we can identify 



mi 



m^ 



2yAt 



m2 



rrib 



V2g„ 



(5.32) 



which leads to the ratio given by 



mt 



m.r, 



W 



iir" 



q2 ■ 



(5.33) 



In the Bogomolny limit, the coupling A is solely given by the Yang-Mills coupling 
constant g as X = Sn'^/g^. The value of A in the Bogomolny limit is smaller than that 
in the London limit, but it is still large A ~ 20 even for g ^ 2, see section |^. 
In the Bogomolny limit mi = m2, the mass satisfies 



m,i = m,2 



2 2 
m = 7 



2r]^ 



(5.34) 



The string tension in the Bogomolny limit is given by 




(5.35) 



since \n{myml) = ln(l + ^1 - 47/77^) - ln(l - ^1 - 4r/V7^) = 2^1 - 4r/V7^- In 
this limit, the monopole condensate reads 




(5.36) 



In the London limit, m,2 reduces to the mass mf, of the dual gauge field and m,i to 
the diverging mass m,^ of the Higgs field. 



2 2/ 2\ 2 '/ 2 \ 

m2->?7(=mj, mi-^ —{=m^^ 00). 
Actually, in the near London case of type II, we have 



(5.37) 
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J'g' - 
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2tt 



p-'Kr,Hn\ 



(5.38) 



which diverges in the naive London limit 7 ^ 00. 
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5.3 Type of dual superconductivity 

In the LEET, the ratio 77^/7^ determines the type of the QCD vacuum as the dual 
superconductor, 

If the ratio rf j^^ is in the range [0, 1/4), the QCD vacuum is the type II dual super- 
conductor. Then a = seems to be excluded. If this ratio is zero, the QCD vacuum 
is the dual superconductor in the London limit. On the other hand, if the ratio ^7^/7^ 
approaches 1/4, the dual superconductor is on the border between type II and type I. 
The Bogomolny limit m\ = 7712 is achieved if the ratio is 77^/7^ = 1/2. The function 
/2(«)//i(«)^ is positive for a > ao and it has a peak 0.21 at a = 2.65, although it 
approaches zero slowly as a increases. The negative factor can come from In -^ in K. 
We can reproduce both types of the dual superconductor by choosing /i. depending 
on the ratio Ma/ jJi- 

Finally, we discuss how the above results change if we consider the effect of V[d]. 

2 

V 

^ "■ ^ ^ ' (5.40) 

p^ T]^ V?"* 7"* ' 

Therefore, the effect is equivalent to the following replacement of the parameters, 

V'^\v^ -,^ -,--,■ (5.41) 

6 7* T]^ 7^ 

Consequently, mi,m2,x are modified as 

ml 2 = , "^^ ,, fl±-J^--V (5.42) 

1 772 

^ = 27F=1- f^'^^) 

For ml 2 to be real and positive, rj^/'y^ must be in the range [3/4, 3). The type II the- 
ory lies in between the London limit rj/'y -^ and the Bogomolny limit rj^/'j^ — > 3/4. 
Therefore, the dual superconductivity of type II is excluded. This result suggests that 
the dual superconductivity of QCD is of type I or on the border between type I and 
type II. Therefore, the correction is important to determine the type of dual super- 
coductivity. The above result will be compared with the lattice results in section 0. 
A simple estimation is as follows. When a = l,cr = 1, rj and 7 are 

27r2 




7^ = ^K'/'Ml = 10F{g)'/'Ml (5.44) 
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and the ratio is 

Fig). (5.45) 



r/^ 27r2 3 ^^ 27r2 3 



74 A^^2 5 ;y^2 5 



l + TT^'-TT^rff' In- ^ 



11 27r2 V 167r^ 



The function F = F{g) is monotonically decreasing in g and goes into the negative 
region for g > 9.6 which is beyond the reach of the analyses of this paper. Finally, the 
string tension in SU(2) Yang- Mills theory is expressed as a function of the Yang- Mills 
coupling constant g as 



where -^((7) := ^^^K{g). More details will be given in a subsequent paper. 

6 Confining string theory and string tension 

Finally, we derive the confining string theory as a LEET of Gluodynamics. The final 
expression of the string action indicates the area law of the Wilson loop or the linear 
static interquark potential where the string tension is represented as a proportional 
constant. It is an old idea that the confining phase of gauge theories can be formulated 
as a string theory, see a review[0]. Especially, the large- A^ QCD might be exactly 



reformulated as a string theory ||45[]. For earlier approaches of QCD (gluon) string in 



the last century, see the references ^, ^\, ^, ^ Q ^ ^, ^, |2 



We can decompose the phase variable ^ as 6* := 6''' + ^* where 9"^ is the regular 



piece and 6^ the singular piece. Then it is shown [Q that the integration measure 
V6 over the field 6 factorizes into the product of measures, i.e., V6 = V6^V6^. The 
singularity of the phase of the scalar field just takes place at the string world sheet. 
The location of the world sheet, x = x(cr), is determined by the condition (j){x{a)) = 
which implies |0(a;)| = on the world sheet where the angle 6{x) is not determined 
uniquely. The singular piece 6^{x) describes a configuration of the vortex string whose 
world-sheet obeys the equation, 

-^e^^'"^(9,9. - d,dp)9'{x) = Q,,{x). (6.1) 

On the other hand, the regular piece 9^ describes a single-valued fluctuation around 
the string configuration just mentioned above. From the correspondence ( p.l[ ), the 
integration over 6'* corresponds to the integration over the world-sheet x^{a) of the 
string, so that the integration measure is transformed as 

Ve'{x)^Vx^{a)J[x], (6.2) 

where J[x] is the Jacobian for the change of the integration variables.]^ 

^-'^The Jacobian J[x] has been evaluated in |Q for the world-sheet of the closed string. The 
Jacobian exactly cancels the conformal anomaly. Such a possibility was already suggested in the 
paper pC| |. Therefore, the vortex string might be self-consistent in _D = 4 space-time dimensions at 
least in the long-distance limit. See also reference pfl, p% p§, |5q |. 
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We return to the expression of the VEV of the Wilson loop ( p.ll8| ) for the action 
m, i.e., 

Sapegt = {h, AD^h) + 2Jgp-'K'/'{h, 6). (6.3) 

Performing the Gaussian integration over the Kalb-Ramond field h^^, we obtain the 
an action written in terms of the vorticity tensor,, 



S,s = J d^x J d'y{2Jgp~'K'/^fQ^,{x) 



X 



X 



A — 7712 A — m( 



(x,y)e^^(i/). 



(6.4) 



for the expectation value of the Wilson loop. 



{W{C))ym = Z,,[C]/Z,M Z,,[C] = jVx^{a)J[x] exp{z5,,[x]}. 



(6.5) 



Here m2, "^25 X ^'^^ the same as those defined in section ^ 




X 



1 1 
m\vn2 

2 2 

m\ — 7712 



1 - V/7^ 



(6.6) 
(6.7) 



Let a = (a , a ) be a two-dimensional coordinate on the world sheet x^ = x^(o"). 



Then the infinitesimal surface element. 



dS^,y{x{a)) = ^Jg{a)t^iy{a)d'^a, 



(6.^ 



is expressed by the determinant g{a) = det ||5'af,(cr)|| calculated from the induced 
metric tensor of the surface defined by 



gab{<j) = daX^{a)dbXf,{a) 



with the derivative. 



da 



d 

da"" 



(a = 1,2), 



and the so-called extrinsic curvature tensor of the surface, 

,a6 



tf,u{<y) 



9[^) 



--daX^{a)dbX^{a) 



(6.9) 
(6.10) 

(6.11) 



Given the theory with the action of the form. 



^=(6, 



K 



A^ — TH? 



0), 



(6.12) 
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the low-energy limit is obtained by performing a derivative expansion of this action. 
The derivative expansion is equivalent to an expansion in powers of 1/m. This proce- 
dure is well know in the literatures, see e.g. [^ ^, |B^, ^, ^ 0, |B^ and Appendix |^. 
The result is the Nambu-Goto action with a rigidity term, 

S,, = a°, j^ (fa^ +^J (fa^g'^'dat^.d^t^" + k^ J Sa^R + • ■ ■ . (6.13) 

The string tension of the Nambu-Goto term is given by 



a 



St 



f^ J. rn 



(6.14) 



where -^'0(2;) is the modified Bessel function and A is the ultraviolet cut-off. More- 
over, it has been shown that the coefficient of the extrinsic curvature term is a negative 
constant, 

11 . ^ 

— = <0, 6.15 

ao 12877 ' ^ ' 

which is independent of k, m, A. 

In the naive confining string theory based on the action ( |6.13|) , the string tension 
diverges 0"^* ^ 00 as A ^ 00, since -R'o(O) = +cxd. This pathology can be automati- 
cally avoided in the confining string theory derived in this paper. It is easy to see that 
the rigidity term cancels in the confining string theory ( |6.4[ ). Therefore, the action 
( |6.4| ) is cast into the confining string action. 



Scs = <yst I d''a^+ — f d^ay^g'^'dat^^dbtf^" + Kt j d''a^R + 
Js ao J J 



(6.16) 



with the string tension, 

ast = {2Jgp-'K 



An 



-«<t)--»(t) 



(6.17) 



Note that the asymptotic behavior of the modified Bessel function Ko{z) for z <^ 1 
is given by 

z 2e~'^^ 
Ko{z) = -i'jE + In -) = In , (6.18) 

2 z 

with 'jE being Euler's constant 'Je = 0.5772 ■ ■ -. Thus, for sufficiently large A, we 
obtain the A- independent finite result Q 



(6.19) 




^^In the setting up of this paper, the role of the ultraviolet cutoff A is played by the Higgs mass 
itih/V^ = mi, as shown below. 

^•^The string tension is a free energy per unit length of the string. It is well know that the free 
energy has a logarithmic dependence in the Ginzburg-Landau theory. The London limit corresponds 
to mi — + cx). Hence, the string tension reduces to the expression, ast = ["^J 9P~^ K^^'^Y ^Kq (^) , 
since it'o(mi/A) -^ 0. 
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where J = 1/2 for the fundamental and J = 1 for the adjoint quark (sources) in 
the case of SU{2), and J = 1/3 for the fundamental quark in the case of SU{3). 
The string tension just obtained agrees with that obtained from the monopole action 
(|5.18|) . The mass mi can be viewed as a dual Higgs mass rriH and l/rrii corresponds 



to a finite thickness of the string. The limit nii — > oo corresponds to the London 
limit, the thin string. This is consistent with a fact that the naive confining string 
theory with the string tension (|6.14 ) is obtained from the DAH model in the London 



limit. 

When the theory is near the London limit mi ^ m2, the expression of the string 
tension reduces to 



>st 



'i^,-^K,nll.(^^)\ (6.20) 

Stt \mb J 



which is similar to the result obtained by Suganuma, Sasaki and Toki |Q. In this 
case, it is again cast into the form, 

<yst = ^^p-'KmlKo (V2^) , (6.21) 

477 V mn/ 

which agrees with the result of Suzuki P] up to a numerical factor. Thus our results 
agree with those predicted based on the hypothetical DGL theory. This fact also 
supports that the DGL theory is one of the LEET's of Gluodynamics. 
The coefficient of the rigidity term is a negative constant, 

an^ = -J^g^p-^K- < 0. (6.22) 

n 

Finally, the k is a positive constant, 

^ = \j'g'p''K^- = -la,' > 0. (6.23) 

String theory was originally developed as dual resonance models to explain hadronic 
physics. It was almost abandoned after the invention of QCD and the discovery of 
asymptotic freedom. Nevertheless, string theory might be useful by offering an alter- 
native but tractable method of solving the strong coupling problem at long distance 
or low energies to which QCD has not yet given the definite answer. Indeed, it is 
known that ordinary strings have the notorious troubles such as tachyons and confor- 
mal anomaly (critical dimensions). However, they disappear asymptotically at large 
distance as shown e.g. by 01esen|]SS|. In other words, a string theory becomes a 
consistent model at large distances, although it is in a strict sense inconsistent in four 
dimensions. Therefore, the large- distance QCD can be described by a string model. 
Of course, such a description breaks down at some distance, since QCD does not con- 
tain tachyons and is Lorentz invariant in spacetime dimensions D less than or equal 
to four. In view of this, it should be worthwhile to mention that the static potential 
for the Nambu-Goto model (i.e., bosonic string model with Nambu-Goto action) has 
been computed to leading order in the 1/D expansion by Alvarez |^ and exactly in 



the whole range of R by Arvis [^ in arbitrary dimension D, 

V{R) = a,,{R^-Riy^\ R, = J^^- (6.24) 
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In the long-distance region R > Re, V{R) is expanded into 

V{R) = astR - ^-^^ + (1/^')- (6.25) 

For the large R, therefore, the static interquark potential is given by the linear poten- 
tial where the string tension agt is the proportional coefficient. Moreover, the static 
potential has an additional long-distance Coulomb term which agrees with the earlier 
observation of Liischer, Symanzik and Weisz |^. The long-distance Coulomb term 



is the universal term depending only on the dimensionality of spacetime due to the 
transverse displacement xt of the string with fixed end points, under the assumption 
that Xt is the only relevant dynamical variable at large distances (This is not the case 
for the superstring |5^). The second term should not be confused with the short- 



distance Coulomb potential which is consistent with the asymptotic freedom of QCD. 
We see that the expression of V(R) loses the meaning at short distance R < Re. 

In this paper we have shown that the Nambo-Goto action can be a piece of the 
effective string action as a LEET of QCD. Therefore, QCD should possess a long- 
distance Coulomb potential. If the string theory which is capable of describing the 
whole energy range of QCD exists, the string theory must reduce in the long-distance 
limit to the above string theory derived in this paper. 

The large N expansion can give another link between QCD and string theory, as 



pointed out by 't Hooft|^. The leading order of this expansion is some kind of free 
string theory that has yet to be identified. In a free string theory, the surfaces in the 
sum should be dominated by smooth surfaces with no surface tension and without 
self-intersections. The Nambu-Goto model describes fundamental string without a 
transverse extension. The string describing color electric fiux tubes in QCD must 
be thick strings with a fundamental transverse length scale mjl (The London limit 
TTiH — *■ cxD corresponds to a thin string) |^. If so, the string action should be 
responsible to the bending rigidity due to the finite width of the string. In order 
to take into account these features, the string model with the extrinsic curvature 



(the so-called rigid string) has been introduced by Polyakov [^ and Kleinert |]48 
The extrinsic curvature stiffness was expected to suppress the crumpled surface with a 
large number of self-intersectoins. However, it turns out that the new term is infrared 
irrelevant, see |^ and references therein for more details. 



Recently, new string theories (so-called the confining string theory) of describing 
the confining phase in gauge theories were proposed by Polyakov [^] and by Kleinert 
and Chervyakov |UT|. The confining string theory has a non-local interaction between 



world-sheet elements and a negative stiffness. The confining string theories are very 
promising, since they seem to solve all the problems of rigid strings. Especially, a 
negative stiffness is crucial in order to match the correct high-temperature behavior 
of large A^ QCD [^2|, |6l[]. The confining string theory can be explicitly derived for 



Abelian gauge theories, compact U(l) gauge theory |Q, Abelian Higgs model[^ 
and so on, see a review |6^. In this paper we have derived the confining string with 
a negative stiffness directly from QCD at least in the low-energy regime. This is 
performed by integrating out the antisymmetric tensor field in the improved version 
of the APEGT (originally proposed by the author in the paper |Tl|) which was derived 
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directly from QCD. The action realizes explicitly the necessary zig-zag invariance of 
confining string |^, |70 . 



7 Parameter fitting for numerical estimation 

For the numerical estimation of physical quantities, we can use the following values 
which seem to be mutually consistent. 

7.1 Dual Ginzburg-Landau theory 

We use the values suggested in [|| For the dimensionful quantities, 

m^ = V2gmV = V2—V ~ 0.88GeV, mn = 2V\v ~ 18GeV, v ~ O.lGeV, (7.1) 
9 

and for the dimensionless coupling constants, 

a, := ^ ~ 0.24 {g ~ 1.7). A ~ 8 x lOl (7.2) 

These values are consistent with the string tension, 

(T,t~(0.42GeV)2~0.18(GeV)2. (7.3) 

The off-diagonal gluon mass obtained by Monte Carlo simulation |2^ for SU{2) is 



Ma = 1.2GeV. (7.4) 

7.2 Monopole action 

In the paper |2^ the following lattice monopole action was adopted, 

S[k]= J2 k,is)D,{s-s')k,is'), (7.5) 

S,s',fl 

where Dq is parameterized by three parameters as 

Do{s - s') = a6s,s' + PAl\s - s') + 7Ai(s - s'), (7.6) 

with the lattice Laplacian Al{s — s') := —dd'. This leads to -Do(p) = « + P/p^ + 7P^. 
Its inverse is 

where 

ml + ml = d/'y, rn\rn\ = (3 / -^ , n \= ^~^ / {m\ — rn\) . (7.8) 

The string tension is obtained as atot = <^d + <^q with 

ad = -n\n — , (7.9) 

2 m2 
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and aq being negligible. The results of ||2^ are 

mi = 1.0 X 10^ m2 = 12, k = 4.83, (7.10) 

and the parameters of the monopole action are obtained as 

a = 0.207(0.435), (3 = 2A9, 7 = 2.07(9.15) x 10-^ (7.11) 

and 

aphys = (0.44GeV)2, aphys/crd ~ (1.4)-^ ^ o.51. (7.12) 

Note that our parameterization is 

ml + ml= 7V^^ mlml = 7^ ml - mj = (7V^^)\/l - V/7^- (7-13) 

Hence the correspondence of the parameters in our theory to the lattice result [E^ 



are given by 

ry2 _ p/^ = 10, 7^ -^ ^/7 = lOl (7.14) 

Hence we obtain 

mi = 10^ ma = 1, x= 10- (7.15) 

7.3 Confining string 

The confining string action has the following parameters (see section 2.1 of the paper 



a ^ 0.2(GeV)^ Kt = 0.003, — ^ -0.005. (7.16) 

tto 

They are obtained from the so-called correlation length of the vacuum and the gluon 

condensate, 

Tg ^ 0.65{GeV)-\T-^ ^ 1.5GeV), ^^((^^fj^) ~ 0.038(GeV)2. (7.17) 

8 Conclusion and discussion 



we 



In this paper, by improving the strategy suggested in the previous paper [^ 
have derived three equivalent LEET's of Gluodynamics, i.e., dual Ginzburg-Landau 
theory, magnetic monopole theory and confining string theory. Although each of 
them has already been proposed and analyzed as a LEET by various authors, we 
have given a first-principle derivation of these theories directly from Gluodynamics, 
i.e., Yang-Mills theory. In other words, we have shown their equivalence in this paper 
as is obvious from the derivation. Especially, we have shown that the monopole 
condensation occurs due to non-zero mass of off-diagonal gluons. 

The very origin of our nontrivial results is reduced to quantum correction to 
the diagonal fields a^ and -B^jy arising from the massive off-diagonal gluons and off- 
diagonal ghosts. This is a novel viewpoint examined in this paper in analyzing the 
low-energy Gluodynamics. In the conventional analytical approaches, the off-diagonal 
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components are completely neglected from the beginning by virtue of the infrared 
Abelian dominance. However, we notice that it is the off-diagonal components that 
convey the characteristic properties of the original non- Abelian gauge theory to the 
LEET written in terms of the diagonal field alone. The reproduction of the (3 function 
of the original Yang-Mills theory in the LEET is a good example of this fact. 

As a mechanism for mass generation of the off-diagonal component, we have used 
the ghost-anti-ghost condensation caused by the quartic ghost interaction in the mod- 
ified MA gauge (section |2!5| ). To author's knowledge, any other analytical derivation 
of off-diagonal gluon mass is not known for the Yang-Mills theory in the MA gauge. 



However, the following steps after section ^TS] can be performed irrespective of the 
origin of the off-diagonal gluon mass, once we regard the massive off-diagonal gluons. 
Therefore, even if the off-diagonal gluons become massive due to other mechanism, 
we obtain the same LEET's as given in this paper. 

The LEET's of the SU(N) Yang-Mills theory have been obtained by way of the 
APEGT. The APEGT is bilinear in the diagonal fields, a^ and B^y. This result 
is regarded as the leading order result of the large N expansion. The (higher-order) 
correction terms are suppressed by a factor A^~^. Since we have obtained the confining 
string theory by rewriting the APEGT, this is consistent with the claim that the 



Yang-Mills theory in the large N limit is equivalent to a certain string theory ||5 
Therefore, our result may shed more light on the relationship between the gauge 
theory and the string theory. However, we have used the derivative (or low-energy) 



expansion or weak field approximation p5| to derive the confining string theory. 



Therefore, it is not clear where the multi-valuedness of the confining string action 
comes from. Multi-valuedness is considered as a reflection of the compactness of the 
residual Abelian gauge group which is embedded in the original compact non- Abelian 
gauge group. This problem will be discussed elsewhere. 

By way of LEET's, we have succeeded to calculate the string tension of QCD 
(gluon) string. The non-zero value of the string tension implies area law of the 
Wilson loop, i.e., quark confinement. However, the string tension obtained in this 
way depends on the parameters p, a which were introduced to rewrite the Yang- Mills 
theory into the dual Abelian gauge theory. It is possible for the string tension to be 
independent of the renormahzation scale [i. In a subsequent paper [^, we will discuss 
in detail the issue how those parameters are determined within the same framework 
as that proposed in this paper. We will give a simple model (toy model for the gauge 
theory) in which the role of parameters introduced in the auxiliary field formalism 
is clarified in more tractable way. We will also give more quantitative argument 
so that the LEET's derived in this paper can reproduce the experimental values of 
physical quantities and predict unobserved quantities, e.g., the glueball mass. We 
hope to discuss the spontaneous chiral symmetry breaking which is expected to occur 
simultaneously with quark confinement. 
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A Useful formulae 

A.l Structure constants 

Note that 

where C2 is called the quadratic Casimir (operator). This implies that 

since the structure constant with two and three diagonal indices are zero, f'^^'^ = = 
p^ according to /^^^ = -2^/^ii{[T^,T^]T^}. Combining the identity, 

with ( [A.ID leads to 

jacdjbcd ^ ^^^ _ 2)5'^^ (A.4) 



A. 2 Differential forms 

For the p-form. 



uj:=-iu^,...^^dx^' A-'-Adx^''^ (A.5) 



the dual form *lj in four dimensional Minkowski space is defined by 

1 
(4-p)! 



*uj: = ,, _ * cu^,...^,,^da:^^ A ■ ■ ■ A rfx^^-^ (A.6) 



The identity, 

fc/xi---Atpai-«4-pt y P-y^ >^/-"[/3i /34-p]' \-^-°) 

leads to 

* *uj = g-\-lYu, (A.9) 

since 

,m-A«p/5i-/34-p ^ g-\^...^^p^...^^_^, g := det{g,,). (A.IO) 

Note that g = —1 ioi Minkowski spacetime with a Lorentz metric, while g = 1 for 
Euclidean space. 
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A. 3 Integration formula by dimensional regularization 

Define 



D 



In Minkowski space, 



i{2n\ 



— ln(m^ + 2p- k- k'^) 



Tie -2) 



2 I 2\2-e 



(An] 



2~e 



(m + p ) 



d^k 



i(27r)^ (m2 + 2p ■ A; - P)" (47r)2-T(a) 



r(e + a - 2) 2 , 2\2-a-e 



d^k 



k,. 



i(27r)^ (m2 + 2p ■ A; - A;2)« (47r)2-T(a) 



'pM +v) , 



(A.ll) 

(A.12) 

(A.13) 
(A. 14) 



/i 



ci^A; 



k^ky 



(27r)^ (m2 + 2p ■ A; - A;2)» 



r(e + a - 2)p^p,K + r )'"""' - r(e + a - 3)-(7^,(m' + p 



2 I „2\3-a-e 



(47r)2-T(a 

where m? is understood as m? — i6,6 > 0, i.e., Im{m'^) < 0. 

A. 4 Gamma function 

The Laurent expansion of the Gamma function is as follows. 



r(e) = --7^ + 0(e), 

r(e-l) = -- + ^E-l + 0{e), 
e 

r(e-2) = i-7i. + ^ + 0(e), 

1 " 1 

T{t-n) = ^-^ --1E + Y.T. 



(A.15) 



n\ 



0(e) 



(A.16) 
(A.17) 
(A.18) 
(A.19) 



where 7^; is Euler's constant 7^; = 0.5772 ■ 



B Derivation of a version of non-Abelian Stokes 
theorem 



In an expression of the non-Abelian Stokes theorem, 

W{C) = J d^MciO exp [tg J^^ dS^^fl 
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fjLV 



x] 



(B.l) 



the argument of the exponential is rewritten as follows. The target space coordinate 
of the surface spanned by the Wilson loop C is denoted by x{a) where a is the world 
sheet coordinate. The surface integral is rewritten in terms of the vorticity tensor as 

^1^ dS^'''{x{a))f^,{x{a)) = j d^xQ>'''{x)f,,{x) (B.2) 

:= (e,/) = (*e,*/) 

= {*Q,/\-\d5 + 5d)*f) 

= {*Q,A-^d6*f) + {*Q,A-^6d*f) 

= {6A-UQ,6*f) + {Q,*A-^6*6f) 

= {6A-UQ6*f) + {Q,A-^d6f) 

= {6A-UQ,k) + {A-^6Q,j), (B.3) 

where k := 6 * f is the magnetic monopole current and j := 6f is the electric current. 
Assuming the absence of the electric source j = 0, we obtain 

W{C) = JdficiO^W [^g{E,k^)+^g{N,J^)] , (B.4) 

where N is one-form defined by 

E := *dQA-^ = 6 * eA-\ N:=6QA-\ (B.5) 

with the components, 

: ]^e^''P^d:jdSQp.{y)A-\y-x) (B.6) 



S^fx) 



1 
-e 



2_ d:j^d'Sp^{x{a))A-\x{a)~x), (B.7) 

N^{x) = d: J d'yQ^%y)A~\y - x) (B.8) 

= -d^ I d^S^"'{x{a))A-\x{a) - x). (B.9) 

C Calculation of the vacuum polarization for ten- 
sor fields 

We shall evaluate the vacuum polarization of the tensor field B. The contribution 
from Fig. ^ is calculated from 

xD,,,,{p + k)5''^'-[-2agr-'-I^p^,,„,l (C.l) 

where 

Ifj.u,aP '■= -{giiaguP — gfifSgua), (C2) 
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and the the off-diagonal massive gluon propagator D"^^{k) = 5°'^D^v{k) is given by 



D^u{k) :-- 



k^-Ml 



k^-M^ 



9tiu -{I -a) 



rCa'^u 



P - aMl 



kaku \ knki, 

9,^u - ^TTT I + 



M2 



M2 A;2 - aM2 



1 

M2 



M^Qf^u - k^K 



+ 



k^ky 



k^-M^ k^-aM^ 

Hence the additional term to the APEGT is given by 

d'^k 



(C.3) 
(C.4) 
(C.5) 



^(1)^ APEGT 



{2n) 



*Bumzuk)*Bu-k). 



Using the identity, 



we obtain 



^fj.iy,al3{k) -^ 



^did2 ricidi^ciC2 pc2d2 _ ricidi pcidi _ q ^ij 



(C.6) 



(C.7) 



-2agf 



Co5'n, 



2 
2a^g^C25'^ 



fiu,piai'iaf^,P2a-2 



' fiu,picri ^ af3,p2CF2 



(2^ 



D^,,,{p)Dp,p2{p + k) (C.8) 



d'^p 



71^4 -pi^,pi(T-i_^fJip,p2(J2 / /OttM 

M^9aia2 - PaiPa2 M'^9piP2 " (P + ^)pi (P + ^)p2 



p^-M^ 



{p + A;)2 - M2 



M'g, 



Pa,Pa2 {P + k)p, {p + k) 



P2 



p2 _ M2 {p + kf - aM^ 



Pa^Pa2 M^9pip2 -{P + k)p, {p + k) 



P2 



p2 - aM2 {p + kf - M2 

PaiP<T2 {p + k)p^ {p + k)p. 



p'^-aM'^ {p + ky-aM^ 



(C.9) 
(CIO) 

(c.ii) 

(C.12) 



Using the Feynman parameter formulas, two denominators are combined into one 
denominator, e.g. 



1 



1 



dx- 



1 



p2 - aM^ {p + ky - /3M2 Jo [p2 + 2xk-p + {xk^ - {xp - xa + a)M2}]2 ' 

(C.13) 
we have 



^pu,al3{k) 



d'^p [M'^9aia2 - PcnPa2] [M'^9piP2 ' (p + k) p^ (p + k) 



P2J 



{2nY 



[p2 + 2a;A;-p + {a;A;2-M2}]2 
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(C.14) 



(27r)4 [p2 + 2a;A; ■ p + {xfc2 - [xa - x + l)M^}Y 

. d^P \Pa1Pa2] [M^9piP2 -{P + k)p, {p + k)p^] 
(27r)4 \p2 + 2xk-p+ {xk^ - {x - xa + a)M^}]^ 

r d^p [Pa^Pa2][iP+k)p^{p + k)p^] 

{2-kY [p^ + 2xk-p+ {xk^ - aM^}y 



(C.15) 
(C.16) 
(C.17) 



The momentum integration can be performed by making use of the formula in Ap- 
pendix ^ where e := 2 — D/2. After straightforward but tedious calculations, we are 
lead to 



^tiu,al3{^) 



X I T{e)[M^ - x(l - x)k^]'' + M-^T{-1 + e)[M^ - x{l - x)k 
-]-M^^T{-l + t)[{a + (1 - a)x}M'^ - x(l - x)k 



,2ll-e 



' 2ll-e 



--M'^T{-1 + e){{a + (1 - a)x}M'^ - x{l - x)k 



.2ll-e 



(to) 



I - M"2r(e)[M2 - x{l - a;)A;2]"^(2a;2 + 22; + 1) 



r 2ll-€ 



,2ll-e 



--M-^r(-l + e)[M2 - x{l - x)k^' 

+M~^T{e)[{a + (1 - a)x}M^ - x(l - x)k^]-'{x'^ + 2x + T 

+ -M-^r(-l + e)[{a + (1 - a)x}M^ - x{\ - x)k 

+M-2r(e)[{a + (1 - a)x}M2 - x(l - x)/^^]-^^' 

+ -M-^r(-l + e)[{a + (1 - a)x}M2 - x(l - x)A; 

-^M-^r(-l + e)[«M2 - x(l - x)A;2]i"4, 



(C.18) 

(C.19) 

(C.20) 
(C.21) 



where the fourth term does not give the term proportional to Ip,v^ap. Here we have 
used the following properties of the projection operator, 

(C.23) 



J^ p,u,al3 
'^ fjLV,pa^ pa,al3 



= I 



pu,a(3i 



-l^pu,pial^ J-ali,p2<T — ~:\f^pi^a9uli f^p,f^p9va l^ul^a9fJ.fi "T kyk^Qpa) 



— -jk {I - P)pu,af3, 



(C.24) 
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where we have introduced 



-i /ii^,Q/3 • ^y-'- ^a-'- u/B -'- fjbli-'-ua} ^ -'- fiv ■ Q^u 



nun 



fini; 



k^ 



Using the Laurent expansion of the Gamma function T{x), we obtain 



Iil...,ik) 



(C.25) 
(C.26) 



-2C2cr'g' ., 






S^-' Ln.nFI^ 



^i^,a/3' 



dx[a + (1 — a)x] 



167r2 



5'n, 



lMU,al3 



dx< 



k^ ■ 



In 



M2 



q; + (1 — Q;)a; — a;(l — x) 



M2 



In 



jo; + f 1 - ajxf^— 



PI 

ft 




— a;(l — a;) 





+ (7ij -ln47r - 1)(1 -a)(l - x) -7E + ln47r 
+ 



X 



2C2cr^g^ ,ii 1 

-2 



167r2 



5^^1A;2(/-P)^,,,^ 



M 





-In 

— M-2 / dx 
2 io 

— M-2 / ^^ 



f dx{2x^ + 2x+l)l In 



M2 



x(l — x) 






|a + (1 - aja;}^ - a;(l - x)^r 



1 — a;(l — x) 



M2 



X In 





'^ / (ix 
Jo 



"-^(1-^)]^ 



In 
In 



^-x(l-x) — 



a- 



q; + (1 — a)a; — x(l — x) 



M2 



P' 
a;(l -a;)^r 

/i2 



{a + (1 - a)a;}^-^^ — x{l-x] ^ 



+ 0(e) 



(C.27) 



M2 ,. .p- 

Here note that the divergent term e~^A;^(l — P) does not exist. In the neighborhood 
of /c^ = 0, i.e, in the low-energy region such that /c^/M^ -C 1, 



n!i„,(A;) 



/ii^,q/3\ 



2^2 



■5'^U, 



fiu,af^ 



2CW_g 
167r2 



1 A;2 



k' 



e '^^ + /o(«) + /l(")]^^ + /2(«J^ 



A;2 
/io(a) + /ii(a)^ 



O 






(C.28) 



where^"^ 
/o(a) 



1+", M2 1 

— - — In — — + (7^; — m47r — lj(l — aj- — 7£; + ln47r 



^^^ 



^* The piece proportional to 1 — P leads to the perimeter law. Therefore, it is neglected when we 
obtain the string tension. 
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- / dx [a + {I — a)x] \n[a + (1 — a)x] 
Jo 

■ln^ + (7i5-ln47r- 1)(1 - a)- - 7^5 + ln47r 



2 
1 



1 — a 



a 



1 1 



— Ina H a^ 

2 4 4 



/i(«) 







(ia; a;(l — x){l + ln[Q; + (1 — a;)a;]} 



(C.29) 

(C.30) 
(C.31) 



1 1 fa^ a^ I ,_, , 

— h 7 7T7 < — In a 1 1 + a 

6 1-a 3| 3 9 9 ^ ^ ^ 



Y^""-T + 4 



+a{a\na — a + l)>, 



/2(«) 



(ia; 



x^(l — x) 
1 



1 3:^(1-0;)^ 
2a + (1 — a)x 



30 d - a) 



1 1 1 3 1412,' 
a H — q; a a ma 

24 3 3 24 2 



(C.32) 

(C.33) 
(C.34) 



and 



ho{a) := / dx< a\na + [a + {1 — a)x]\n[a + {I — a)x] 

Jo 2 



hi(a) 



(2x2 + 2x + 1) ln[a + (1 - a)a;] ^ , 

^ , f / 9 N (1 — «)a;(l — a;)2 

l^ a + (1 — a)x 

H — x(l — x) In a — x(l — x) ln[a + (1 — a)x] >. 



(C.35) 



(C.36) 



For a = 1, the results are greatly simplified; hi = and /o = — In ^— 7£; + ln47r, /i 
1/6, /2 = 1/60. 



D Manifest covariant quantization of the second 
rank antisymmetric tensor gauge field 

D.l Second- rank antisymmetric tensor gauge theory 

We discuss the gauge fixing of a second-rank antisymmetric tensor gauge field A^^, 
whose Lagrangian is given by 



1 



u Apa\2 



o 
This Lagrangian is invariant under the hypergauge transformation, 



(D.l) 
(D.2) 
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In order to fix the gauge, we adopt the gauge fixing condition for A 



liui 



d^A 



fJ,U 



0. 



:d.3) 



Then the gauge fixing (GF) and Faddeev-Popov (FP) ghost term is obtained based 
on the prescription of Kugo and Uehara [^ 



as 



Ci 



-i6f 



a\ 



C\d^A^,^^B, 



:d.4) 



where we have defined the nilpotent BRST transformation, 



^BC^,[x) -- 

6Bd{x) = 

5bC^{x) -- 

5BBf,{x) -- 

Hence, the explicit form of Ci reads 



idf,d{x), 

0, 

iBf,{x), 

0. 



£i = B'-d^A,, + iCd^ld^C, - d,C,] + -j{B^) 



(D.5) 



(D.6) 



However, the Lagrangian Ci and hence Cq + Ci is still invariant under the trans- 
formation of the vector ghosts C^ and (7^, i.e., 5Cfj,{x) = idfj,9{x),6C^{x) = id^ip{x). 
Note that C^ and (7^ are independent fields and that C^ = C, C"!" = C. We consider 



the gauge fixing conditions, d^C^ 
additional GF+FP term. 



and d^Cn = 0. Therefore, we must add an 



C, 



45} 



d{d^C^ + a2P) 



iSf 



N (d^C. + a^B'-^^ 



(D.7) 



where the nilpotent BRST transformation of the additional fields is supplemented as 

8bN{x) = P(x), 

5bP{x) = 0, 

5Bd{x) = B^^\x), 

6bB^'\x) = 0, 



The explicit form of £2 reads 
£2 = -tB^^^d'^C. 



la^B^^'^P + dd^'d^d - iPd^'C^ + Nd^'B^, 



(D.8) 



(D.9) 



where we have defined a^ := 0:2 — 03. Note that P and 5*^^^ anti-commute. For 
the assignment of the ghost number of each field, see Table. 1. Two vector fields 
CiJLi Cfj_ are two primary ghosts, and three scalar fields d, rf, A^ are three secondary 
ghosts. Three fields -B^, P, B^^' are the Lagrange multiplier fields for the condition, 
d'^A^u = 0, d^'C^ = 0, d^'Cf, = 0, respectively. Thus we obtain the GF+FP term for 
the Lagrangian Cq, i.e., Cgf+fp = -^i + £^2- 
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field rank ghost number 



A 


2 





C 


1 


1 


d 





2 


c 


1 


-1 


B 


1 





N 








P 





1 


d 





-2 


S« 





-1 


A 


1 





C 





1 


C 





-1 


B' 









Table 1: The ghost number of the field with the indicated rank. 
Now all the gauge freedom is fixed. Thus the full Lagrangian density is given by 

^tot = -^0 + ^GF+FP, 

-iB'^^^d^C^, - la^B^^^P + dd^'d^d - iPd^'C^ + Nd^'B^. (D.IO) 

The massless antisymmetric tensor field stands for the massless spin-0 field as a 
physical mode. It is possible to show that all the unphysical modes decouple leaving 
correctly one physical mode |^ . 



The above result is the summary of the results obtained by Townsend |38], Hata 



Kugo and Ohta ^^ and Kimura [^]. The same result can be obtained within the 



framework of the extended theory for the constrained system based on the canoni- 
cal Hamiltonian formalism on the extended phase space, the so-called the Batalin- 
Fradkin-Vilkovisky (BFV)formalism, see e.g. the original papers and a review |7^ . 



D.2 Inclusion of mass term 

Next, we consider the theory of an antisymmetric tensor field with the mass term, 

C^[A] = -\{e,.,.d^ An' - \m\A,;)\ (D.ll) 

This Lagrangian with the mass term is no longer invariant under the hypergauge 
transformation of A^^. However, the invariance is recovered by introducing an addi- 
tional vector field A^ in such a way that 

£^'[A, A] = -he^^^^d- AP'^f - l(mV + d^^u - d^.f. (D.12) 

o 4 
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Actually, this Lagrangian is invariant under the combined transformation, 

6A^^{x) = d^^^{x) - d^^^{x), 6A^{x) = -m^^(x). (D.13) 

Moreover, it has an additional invariance under the transformation, 

SA^,ix) = 0, 5A^ix) = d^coix). (D.14) 

Therefore, we define the BRST transformation of the field A as 

6bA^{x) = -mC^ix) + d^C'ix), (D.15) 

by introducing additional ghost C. The BRST transformation of C and C' is defined 

by 



SsCix) 
SBC'ix) 
5bB\x) 



imd{x), 

iB'{x), 

0. 



(D.16) 
(D.17) 
(D.18) 



The first BRST transformation is determined by nilpotency of (p.l5|) . Other two 
transformations are also suggested from nilpotency of BRST transformation. Thus 
the nilpotent BRST transformation is determined for all the fields. In order to fix the 
gauge freedom, we must obtain the GF+PF term Cgf+fp- For C[^^ := C^' + Cgf+fp, 
the generation functional of the theory is given by 

Z := fvA^^VA^VB^VC^VCf.VdVdVNVPVB^^'^VC'VC'VB'exp i f d^^xC^t . 

(D.19) 
A good choice isQ 



-GF+FP 



46 



B 



C d^'A 



fif 



d,N - aK + Y^, 



+C' 9^ A, 



a 



B' 



d (d^C^ + bC + a2P) 

(D.20) 



B" d^A 



«! 



/lU 



d,N - aA, + -^S 



-iB^^^d^'C^ + bC + a2P) + dd^'d^d + bmdd 



a 



+B\d^A^ + -B') + iC'd^'id^C' - mC^) 



(D.21) 



where ai,a2,a' are gauge fixing parameters and a, 6 are parameters specified later. 
From the naive viewpoint, this corresponds to the gauge fixing condition, d^A^j^ = 0, 
d^C^ = = d^^C^ and 9^A^ = 0. If we integrate over B and B', the sector containing 
A and A reads 

(D.22) 



Z := f VA^^VA^VN exp i f d'^xCl^ 



25 



The author would Uke to thank Atsushi Nakamura 173] for helpful discussions on this Appendix. 
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where 



C = C^[A]--mA^,{d^K-d,K^)--m\d^K-d,k,f 
la l(y.\ 



(D.23) 



since other fields decouple from the relevant sector. Then we perform the integration 
over iV and obtain 

Z = |pVexp|zy'd^x(/:™[A] + ^(9'^V)')} 

X fvA^expfi fd^xC'^[A,A]\, (D.24) 



where 



CT[A,A] 



--mA^^id^A^ - d^^) - -w?{d^A^ - d^ 



4 

+ ^{A^f + ^{d^A^,)A, + ^{d^A^)^-\d^A,). (D.25) 
2ai ai f^ f^ 2ai 



If we choose a = mai, the cross term {d'^A^^)A^ cancels with mA^^{d^Ai, — d^A^). 
Hence, £^[A,A] becomes independent of A field. Thus the A field decouples from 
the theory of A, 



£-[A] = --m\d,A, - d,A,f + ^{A,f + ^{d''A,)/\-\d'A, 
4 lai Zai 

By taking the Landau gauge ai = 0, we recover the partition function, 
Z = J VA^Jid'A^,) exp |z J d*xC^[A] 
Moreover, it is possible to consider a simpler gauge | 76| , 



CA^, + dC 



^GF+FP — —i^B 

= B^A^ + iC^id^C - mC^) - iB'-^^C + mdd. 



(D.26) 



(D.27) 



(D.28) 
(D.29) 



This corresponds to the gauge fixing condition, A^ = and C = 0. The integration 
over N, P and B' is trivial, since C'^^^ does not include them. The B integration leads 
to the constraint 5(A^), 



Z = f VA^^VA^6{A^) f VC^VCf,VdVdVB^^^VC'VC' exp i f d^xCl^ 



where 



Clt = C^'[A^., A^] + iC^{d,C' - mC,) - iB^^^C + mdd. 



(D.30) 



(D.31) 
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(a) ^.'-.i^--.^ (b) (0) /-T,^ 



\ / 



^2 i!L_.'P2 ^2_1^^ ^L_.^2 "^2 L__.^: 

(d) 



'2 




Figure 5: Self-energy diagrams for the would-be Nambu-Goldstone particle ip2- 



When we consider the sector of A and A, the sector described by other fields decouples 
and we obtain 



Z = I VA^,VA^S{A^) exp jz J d'xC^'[A^,, K)} = J ^ V exp [t J d''xC^[A^A 

(D.32) 
Therefore, we recover the original theory which is written by the A field only with 
the Lagrangian C^. Note that the massive antisymmetric tensor gauge theory stands 
for the massive spin-1 theory. 

E Renormalization of dual Abelian Higgs model 

E.l Goldstone bosons remain massless in higher orders 

We examine the self-energy diagrams for the would-be NG boson which are given in 
Fig. i 

E.(0) := (-2,A„)3^(-6,:A.)i/|^^^ = -3A/(™'), (E.l) 
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(a) 



9- 
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(b) 



^1 



(0 



/' ^1 ^\ 



/ ^ u 




^^ ^2 y 



n 



Figure 6: Vacuum polarization diagrams for the 6^ particle. 



S.(0) 



'-2iXv) 



-m^ 



-2iXv)-^ 



d'^k 



2 J (27r)4p-^M2 



-\I{iM'^ 



S,(0) := {-2i\vf j 



(E.3) 
(E.4) 



(27r)4 A;2 - ^M^ A;2 - ni^ 



4A2 



2A - ig^ 



[I{m')-I{^M% 



S^(0) := (-2iAt 



-m^ 



i^g^M) 



d^k 



(27r)4 P - ^M2 



e^'/(eM^ 



(E.5) 
(E.6) 



When ^ = 0, it is easy to check that the total sum of the above self-energy parts 
vanishes. This fact ensures that the Nambu-Goldstone particle remains massless even 
if we include the higher order corrections. 

E.2 Vacuum polarization 

Next, we consider the vacuum polarization for the b^ field. The relevant diagrams 
to one loop are enumerated in Fig. ^. The vacuum polarization tensor for the dual 
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vector field b^ is obtained by summing up the following contributions, 



(27r)4 F - m^ 



:=2fi^ = 2Xv^, 



n(fj(p) := 2tglg^,J 

^P) ■■= {2^9lvgJ{2^glvg^^) J ^^^^^ ^^ _ ^^^ _ ^^ ^^ ^^^^ 



(E.7) 



X 



A; k 






n(rj(p) := 2tglvg^,tg^M J 



(27r)4 P - ^M2 



(27r)4 F - ^M2 



(E.8) 
(2A;-p)^(-2A; + p),, (E.9) 

(E.IO) 
(E.ll) 



-m 



First, the dimensional regularization yields 

d^k I .^ o f d^k 



nJ^Hp) := 



■^9m9lJ.v 



[ d k 1 _ •<. 2 [ 



i(27r)^ -P + 



m" 



'^■^9m9fiv 



Tit -I] 



(4^; 



,2-e 



m 



2\l-e 



(E.12) 
(E.13) 



i2gl,g^,—-^ (-- + 7s - l) "^'(1 - elnm2)(l + eln47r) (E.14) 



i^g^,um'(-K-l + \nm^ 



where we have defined e := 2 — D /2 and 



(E.15) 



N^ := - + ln47r-7ij. 
e 



(E.16) 



and used the Laurent expansion of the Gamma function r(e — 1) 
Second, we want to calculate 



n!S(p) 



:= i^^9lvgJ{2^glvg.^) J ^^^—^^—-^^^ 



kpkfj 



(27r)4-(p-A;)2 + m2 



9tJ.u M^ 



~M2 



A;2 + M2 -A;2 + {M^ 



(E.17) 



The Feynman parameter formula can combine two denominators into a common de- 
nominator. Then the dimensional regularization leads to 

nS(p) 
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^"^ Jo [J i{27r)D [-A;2 + 2xp ■ k - xp^ + xm^ + (1 - x)M^]^ 



d^k 



~M2" 



i{27r)^ [-A;2 + 2xp-k- xp^ + xm^ + (1 - x)^M^]^ j 
^"^ ' dxl r{e)g^4{x^ - x)p^ + xm" + (1 - x)M^]-' 



(E.l^ 



(47r)2-T(2) Jo 

-r(e)a:'^[(a;2 - a;)/ + xm" + (1 - a;)M2] 



2l-e 



£ fiV^r/^2 

'2M2' 



+r(e - l)l^[(a;2 - x)/ + xm' + (1 - a:)MY"' 



2i-e 



+r(e)a:2^[(a;2 - x)p' + a^m^ + (1 - a:)eM2] 
-r(e - 1)^|^[(:^' - x)p' + xm^ + (1 - x)eMT-4. 
Hence, using the Laurent expansion of the Gamma function r(e), we have 

nS(p) 

g^lu {Ne - ln[(a;^ - x)p^ + xm^ + (1 - a;)M^]} 



(E.19) 



(47r)2r(2) io 
2M2 



(ix 



+-7^[(^' - x)p' + xm^ + (1 - x)M2] 



i-N^ - 1 + ln[(a:^ - a;)p^ + xra^ + (1 - a:)M^ 
2M2 



-I 



X 



{iV, + 1 - ln[(x^ - a;)p^ + xrn^ + (1 - a;)^M2]} 



+X' 



2PMPi^ 
M2 



+ ln[(x^ - x)p^ + xm^ + (1 - a;)M^ 



- ln[(a;^ - x)p^ + xm^ + (1 - a;)^M^] 

Using the low-energy expansion, 

ln[(x^ - x)p^ + xm^ + (1 - a;)^M^] 
= ln[xm^ + (1 - a;)^M^ 



0( 



e . 



(E.20) 






xm2 + (1 - x)^M2 



p^+ofe . (E.2i; 



we obtain 

n2)(p) 



(4vr)2r(2) 
/o 






(7^, / rfx <^ AT, - ln[xm2 + (1 - x)M2] + - 



p 



+ \^tdx[xm' + {l-x)M'' 
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X — X 



X - iV, - 1 + \n\xm' + (1 - x)M^] + ^— —-p" 

y xm^ + [l — x)M'^ 

+ -T7^P'^ t dxix"^ -x)l -N,-l + Inlxm^ + (1 - x)M^] I 
2 M'^ Jo [ J 



X J iV, + 1 - ln[xm^ + (1 - a;)^M2 






P 



xm2 + (1 - x)^M2 

+ -^p' /' dx{x^ - x)\n, + 1 - \n[xm^ + (1 - a;)eM2] 1 
2 M'^ Jo { J 



PmPi' 



/" dxx^l + ln[xm'^ + (1 - a;)M^] - ln[xm^ + (1 - a;)^M^ 



M2 JO 

+Oip'). 
In particular, when ^ = 0, we obtain 



(E.22) 



n2](p) 



^2 \ 



3, , m^lnm^ -M^lnM^ 1 m^ 

i(^ + i) ^^zrw^ 8mJ'-i + 2'""' 



-0 



M2 



Third, we can proceed in the similar way to (b) and obtain 

n!rJ(p) 



g^^{N, + l)[--p^ + m^-^M' 



(E.23) 



(E.24) 



(4vr) 



-2^, 



/J,l/ ^ 



dx[xm'^ + (1 - a;)^M^] ln[a;m^ + (1 - x)^M'^] 



+p'^ / (ia;(x^ — x) + p"^ / (ix(x^ — x) ln[xm^ + (1 ^ 2;),^M^] 



1 /■! 

+ -N^p^Pi, — p^Pu / (ia;(4a;^ — 4a; + 1) ln[a;m^ + (1 — x)^M'^ 

O JO 



0(ptJ5.25) 



In particular, when ^ = 0, we obtain 



nS(p) 



^T^I^M. 



(47r)2 



p 



2 1 N, 13 ^, ^\ 2/Ar 3 , 2 

' —buim]+m\]S^^ mm 
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mm 

3 9 3 



0{p'). 



(E.26) 



The remaining quantities, H^"^) and II^'^j are calculated in the same manner as n^"\ 
However, they vanish for ,^ = in the dimensional regularization. 
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Therefore, the total sum of the vacuum polarization for ^ = is given by 



(47r)2 
(47r)2 
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m^- M^ 
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(4vr) 



-m 5f^^ A/", + - - In — - 



Hm 



(4vr) 



P ^M^ 



iV, 13 
3 36 
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+ 6 In 
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+ PmP^ 



-1 +21n. 



m 



N^ 4 1, m^ 

^ H In — - 
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p^ 
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(E.27) 



where m^ = 2\v and M^ = (74"^. 

Finally, the renormalization factors are obtained as 



Zt(MtY - (Mi 



''^' = -(l|5'"'H-l + '"- 



4q2 



(4vr) 



K^m^) 



-(iV, + l) 



m^ln^-M^ln^ 1 m^ 



m2-M2 



M2 



-1 + 2 In 



m 



(47r 



3 

2 



m 



™ m2(iV, + ^-ln^)+0((?^ 



/i- 



(E.28 



F Calculation of the Wilson loop 

In the following we show that the area law decay of the Wilson loop is obtained from 
the result, 

{W{C))yM = exp [-^{2Jgp''K'/y{%, Z^^^S^)} , (F.l) 

where S is the one-form defined by 



with the component. 



E'^ix) 



E:=*deA-^ = 6*eA-\ 



e^-^W: d%Qp,{y)A~\y-x) 



X) 



^e>'^^W:Jj'S^,{x')A~\x' 



(F.2) 

(F.3) 
(F.4) 
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Then the argument of the exponential is cast into the following form, 

= (6, A-iD„i(A)e) - (6, A-'d5A-'D-,'e) 

= (0, A-iZ}^i(A)0) - (50, A-2D^i50). (F.5) 

For the rectangular loop with side lengths T and R in the xi — x^ plane, we take 

Qo.p{z) = 6ai6^4S{z2)6{zs)e{zi)e{R - zi)e{z^)e{T - z,). (F.6) 

Then the Fourier transformation is given by 

Q^p{p) = J d'ze^p{z)e-'P-' 

= ^aiV / dz^e-'P'"' f dz^e-'P*'^ 
JO Jo 

2 _j£i« . PiR 2 _^E4I. . PaT 



^ai^pi—e^ 2 sm- 

Pi 2 P4 



2 sm ■ 



In the momentum representation, we have 

(0,A-iZ}^i0)=|^0„^(p)0«^(-p)[A-iZ},;i](p). 
If we use the formula following |2^ , 



lim 

R^oo 



sin aR 



T:R5{a) 



(F.7) 



(F.8) 



(F.9) 



for large R and large T, then we obtain 

d'^p 



(0,A-^D„,i0) = j j^,{2^fTR6{p,)6{p,)[A-^D;^\p) 



TR 



-TR 



d^p 
(2^ 
d^p 



[A-iD^i](0,p2,P3,0) 
1 



(27r) 



pI + pI + ml P2+pI + ml 



(F.IO) 



Here the logarithmic divergence of the integral is removed by introducing the ultra- 
violet cutoff A as 

d'^p 



(0,A-^D-^0) 



'TR , , 

(27r)2 

/.A2 

-TR lim 



1 



1 



pi+pi + mf pl+pi + mi 



d\p\^ 



A^oo 7o 47r 



-K 



IpP + mf IpP + ml 



_ K , , m? A^ + m^ 
Ti?— lim In-— ^ -^ 

4:71 A^OD A^ + mf 1712 



K , m^ 



TR— In 



4:71 ml' 



(F.ll) 
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So we obtain the A-independent result. The last term in ( [l''.5| ) gives a perimeter decay 
part, since 60 is the boundary of the rectangular surface. The dominant term in the 
large loop is given by the contribution ( [bMlD which exhibits the area law. Thus we 
arrive at the result ( |5.26|) . 



G Derivative expansion of the string 

In this section, we begins with the expression. 



{W{C))ym 



exp 



Sc J Sc 



where 



G^u,p.{x,y) := -2.Pg''p-^K{V[d,]hl,{x)V[dy]hlM)APEGT. 
In this paper, we have obtained the result, 

X 



-2J2^V"'^V,P- 



X 

A — ml A — ml 



We define the Euclidean propagator: 

GUx) = {-AE + m^)-\x,0)= J 
It is written in the form, 



d ^ ik-x ^ 

(27r)4 P + m2' 



Gm.(x) = —^^-Ki{m\x\) = —^m'^ 



47r2 



x\ 



47r2 



ki/e 



(G.l) 

(G.2) 

(G.3) 
(G.4) 



(G.5) 



(G.6) 



where Ki{z) is the modified Bessel function and we have defined the correlation length 
^ by ^ = m^^. (|G.6|) is obtained as follows. Substituting the identity. 



^2 _|_ ^ jQ 



dse-'^'^'^'^'K 



(G.7) 



into ( |G.5|) and performing the Gaussian integration over the four momentum k, we 
obtain 



\Jfrn \ "t^ / 



dse 



.2 / d^k 



-sk +ik-x 



{2ny 



/ dse "*"■ exp 
Jo 

/ ds^r exp 
JO s^ 



X 



(27r)- vv s 



1 r^ . 1 



—sm, — 



X 

4i 



(G.8) 
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The above result ( |G.6| ) is immediately obtained by applying the integration formula 
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r°° 1 
/ dxx'^ exp 
Jo 



7X 

X 



2{(3r^y'^KM(3^) {^(3 > 0, 3fi7 > 0)), (G.9) 



to the case z/ = — 1,/3 = x^/4,7 = rn?, since K^y{z) = K^{z). 
In Euchdean space, 



^ fiu,pa\Xy X J 



2.fg^p-^KI^,^p,K [G^, (x - x') - G^, (x - x')] (G.IO) 

(G.ll) 

(G.12) 



2V,p.[Fi((x - x')') - F2((x - x')')] 



7-1// M2\ t2 2 -2 r^ X 1TiiKi[mi\X X 

Fi{{x-x} ) := -J g p K- 
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'^^2^ 



J, := 2 / rf5'^'^(x(a)) / dSP''{x{a'))I,,^,,F,{{x{a) - x{a')) 

J Sc J Sc 



It is shown [|65| , |6q1 that the derivative expansion in powers of 

C'^:=(a'-a)7e,, 6 := m'S 
leads to 



4e/M^--efM^^"''(9,V)(9,V) 



J. = / t/V^ 

with the moment, 

Mi := I d^ziz^Fiiz^), z"" := (^'^V 



o(ef), 



(G.13) 



(G.14) 



(G.15) 



(G.16) 



where we have used the conformal gauge for the induced metric, gab{cr) = y g{<^)^abi 

hence C"'C''9ab = g~^^^gabz"'z'^ = z'^z'^^ab '■= z^. Thus, the confining string theory 
derived in this paper is characterized by the parameters. 



a = 4 f d^z[mfFi{z^) -m-^F2{z% 
«o"' = -^ldhz^[m^'Fiiz^)-m^^F2{z% 
n = ^j d^zz^\m:[^F^{z^)-m2^F2{z)]. 
By substituting ( |G.12D into ( |G.17D , we obtain 
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X 
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2iT\z\d\z 
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(G.17) 
(G.18) 
(G.19) 

(G.20) 
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where we have used Ki{x) = —Ko{x), Ko{oo) = and introduced the ultraviolet 



cutoff A (-ft'o(O) = oo). Incidentally, the asymptotics of Kp{z) for z > 0, 



K„ 



e-^[l + 0{z-% 



(G.22) 



means that Kp{z) decreases exponentially for large z. After removing the cutoff A, 
the above expression reduces to a finite value, 



_ (2J^ _, fm. 

Air \m2 

The coefficient of the rigidity term is calculated as 

\z\^K,{\z\) foo 



(G.23) 
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J'g'p-'K 



X 
47r2 
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An An 

2 2 

mf 1712 



= -jVp-^K-<o, 

n 
where we have used the integration formula. 



KAjJijtXj -t V ij 



ax] 



2-^a-r (^^^\ r (^^±^ 



{^p > ^u) 



(G.24) 



(G.25) 



Here note that the integral in ( p.24|) is finite and we don't have to introduce the 
cutoff. Finally, the k is calculated as 



6 n 



(G.26) 
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